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FOREWORD 


Opinions,  interpretations,  conclusions  and  recommendations  are 
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material . 

_  Citations  of  commercial  organizations  and  trade  names  in  this 

report  do  not  constitute  an  official  Department  of  Army 
endorsement  or  approval  of  the  products  or  services  of  these 
organizations . 
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Research  Council  (NIH  Publication  No.  86-23,  Revised  1985) . 
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INTRODUCTION 


Subject:  This  project  optimizes  and  evaluates  a  computational  clustering  technique  for 
aiding  the  diagnostician  to  discriminate  malignant  breast  lesions  during  a  post  contrast 
agent  injection  time  series  of  three-dimensional  MRI  images.  Purpose:  Currently,  few 
diagnosticians  are  skilled  in  the  art  of  identifying  malignancy  from  these  images.  It  is 
possible  that  the  contrast  agent  may  diffuse  into  and  out  of  malignant  tumors  with  a 
unique  time  course  signature  discemable  via  a  fuzzy  clustering  algorithm.  Scope:  The 
project  employs  an  existing  fuzzy  clustering  algorithm  previously  used  to  discriminate 
brain  tissues.  Several  computational  parameters  must  be  optimized  for  he  clustering 
algorithm  to  be  useful  in  the  setting  of  breast  lesions,  and  the  usefulness  of  the  technique 
in  discriminating  malignant  vs.  benign  lesions  must  ultimately  be  evaluated. 


BODY 


% 


Background  and  Introduction 

The  Fuzzy  Clustering  Algorithm  Implimentation 

The  fuzzy  clustering  analysis  (FCA)  properties  of  robustness  and  being  unbiased  are 
exactly  the  characteristics  which  make  FCA  an  attractive  method  for  analysis  of  dynamic 
gadolinium  time  course  curves. 

An  advantage  of  using  fuzzy  clustering  is  its  ability  to  separate  the  gadolinium 
uptake  data  into  K  distinct  clusters.  As  we  have  implemented  it,  associated  with  each  of 
the  K  clusters  is  a  temporal  cluster  centroid  and  a  cluster  membership  map.  The  temporal 
centroid  displays  the  temporal  course  of  the  gadolinium  enhancement.  It  is  calculated  as 
the  weighted  mean  of  the  normalized  temporal  data,  where  the  weight  of  each  temporal 
course  is  proportional  to  the  corresponding  membership  to  the  cluster.  Thus  temporal 
courses  that  have  high  membership  to  a  cluster  contribute  strongly  to  the  cluster  centroid, 
whereas  temporal  courses  with  low  membership  do  not. 

The  cluster  membership  map  gives  the  “degree  of  belonging”  of  each  voxel  to  the 
cluster;  it  displays  where  in  the  breast  the  gadolinium  uptake  curve  depicted  by  the 
corresponding  cluster  centroid  occurred.  As  opposed  to  hard  clustering  algorithms,  for 
which  the  weights  are  either  one  or  zero  (i.e.,  belonging,  or  not  belonging),  the  fuzzy 
clustering  algorithm  associates  a  fuzzy  membership  (a  value  between  zero  and  one)  for 
each  spatial  location  (i.e.,  each  voxel)  to  each  of  the  K  clusters.  The  fuzzy  membership  is 
analogous,  but  not  equivalent,  to  the  probability  of  belonging  to  the  cluster.  This  more 
realistically  represents  the  dynamic  processes  that  occur  in  perfusion  breast  data. 

Statement  of  Work  (First  18  Months) 

Task  1.  To  segment  breast  tissue  on  the  basis  of  MR  image  contrast  enhancement  kinetics.  Determine  the 
fuzzy  index  and  threshold  for  membership  in  a  cluster,  (months  1  to  8).  (Status:  complete.  See  methods 
and  materials  computation  seciton.)  Determine  the  “normalization  or  preparation  of  the  time  course  data 
needed  for  presenting  the  images  to  the  fuzzy  clustering  algorithm,  (months  1  to  12).  (Status:  in  progress) 
Using  the  parameters  determined  above,  perform  fuzzy  clustering  on  300  image  data  sets  and  determine  the 
sensitivity  of  the  fuzzy  clustering  for  detecting  breast  abnormalities,  (months  12  to  24).  (Status:  completed 
43  of  the  300  sets.) 


Methods  and  Materials 

Patient  Population 

To  date  we  have  performed  fuzzy  clustering  analysis  on  43  breast  MRI  cases.  Breast 
images  from  seven  women  are  used  here.  All  subjects  were  selected  for  MRI  because 
they  had  a  suspicious  finding  on  mammography  or  a  palpable  lesion  and  were  also 
scheduled  for  excisional  biopsy. 

MR  Image  Acquisition 

Magnetic  resonance  images  (1.5T)  of  the  breast  were  obtained  using  3D  volume  fast 
rf-spoiled  gradient  echo  imaging  scheme  with  fat  suppression.  To  obtain  3D  images  in 
90  seconds,  a  partial  data  acquisition  scheme  with  32  kHz  acquisition  bandwidth  was 
used  which  allowed  TR  times  of  9.9  msec.  The  TE  time  was  2.2  msec  and  provides 
fat/water  opposed-phase  images.  The  flip  angle  of  the  excitation  pulse  was  45  degrees. 
The  fat  suppression  pulse  consists  of  a  frequency  selective  inversion  pulse  applied  at  the 
largest  lipid  resonance  (220  Hz  from  water)  once  every  16  TRs.  The  inversion  delay  for 
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best  fat  suppression  was  optimized  empirically.  The  image  acquisition  consists  of  a  data 
matrix  of  512  x  256  x  32  resulting  in  an  image  with  512  x  512  x  28  matrix.  The  field  of 
view  of  the  images  is  18  cm  with  a  slice  thickness  that  ranged  from  0.7  to  3  mm 
depending  on  the  size  of  the  breast.  The  slices  are  in  the  sagittal  plane.  To  attain 
sufficient  signal-to-noise  for  these  high  spatial  resolution  images,  a  four  coil  receive  array 
was  used.  The  breast  was  lightly  compressed  between  two  plates  which  contain  two  coils 
each.  The  compression  helped  the  breast  to  fill  the  1 8  cm  square  FOV  while  narrowing 
the  breast  in  the  slice  direction.  For  all  studies,  one  pre-  and  four  post-gadolinium  images 
were  obtained  at  90  second  time  resolution. 

Computations 

Regions  of  significant  gadolinium  uptake  can  be  defined  from  the  cluster 
membership  map  of  each  cluster  by  thresholding  the  map  at  a  selected  level  of 
membership.  This  thresholding  process  creates  partitions  (i.e.,  groups)  from  the  clusters. 
A  threshold  of  0.80  to  0.90  typically  performed  well.  Higher  values  would  be  ineffective 
because  the  presence  of  noise. 

Intrinsic  to  the  fuzzy  clustering  algorithm  is  the  fuzzy  index,  m  (m  >  1 .0),  which 
controls  the  "fuzziness"  of  the  clusters.  A  value  near  one  results  in  "hard"  clusters, 
whereas  larger  values  result  in  successively  more  "fuzzy"  clusters.  A  compromise 
between  hard  and  fuzzy  clusters  was  found  to  be  m  =  1.2,  which  has  been  used 
successfully.  This  value  has  also  been  supported  by  studies  of  both  simulated  and  actual 
fMRI  data. 

Results  and  Discussion 

Data  Summary 

Space  limitation  preclude  the  possibility  of  conveying  images  for  the  43  cases 
acuired  to  date.  Typical  results  are  shown  for  cases  representing  a  benign  lesion,  ductal 
carcinoma,  fibroadenoma,  and  lobular  carcinoma. 

Benign  Lesion 

Figure  1  shows  the  results  of  a  first  pass  FCA  in  a  breast  with  a  benign  lesion.  The 
two  clusters  containing  uptake  curves  with  significant  enhancement  gives  the 
parenchyma  and  the  lesion  (clusters  1  and  2).  Clustering  on  shape  of  the  curve 
redistributes  the  pixels  in  clusters  1  and  2  from  the  first  pass  into  four  clusters  (Fig.  2). 
The  pixels  within  the  lesion  are  primarily  in  cluster  1 .  The  centroid  for  cluster  1  shows  a 
peak  in  enhancement  at  90  milliseconds. 

Figure  1 _  _  _ 


Ductal  Carcinoma 

Figure  3  shows  the  centroids  and  membership  map  of  cluster  corresponding  to 
centroid  4  for  three  slices  of  this  breast.  There  are  two  lesions,  one  in  the  first  slice  and 
another  in  the  third.  Partial  voluming  of  both  lesions  in  slice  2  can  be  observed  as  well  as 
some  pixels  that  bridge  the  region  between  them.  The  histopathology  indicated  two  foci 
of  invasive  carcinoma  with  ductal  carcinoma  connecting  the  two  lesions.  The  two  lesions 
are  easily  picked  up  by  the  FCA  as  well  as  the  bridging  DCIS  (ductal  carcinoma  in  situ) 
in  the  intermediate  slice.  A  third  pass  through  the  FCA  on  the  lesion  cluster  gave  three 
distinct  clusters  (see  Fig  4).  One  cluster,  2,  consists  of  the  lesion  while  cluster  1  is 
predominantly  found  along  the  rim  of  the  lesion.  Careful  examination  of  the  borders  of 
the  lesion  under  the  microscope  revealed  that  it  was  bordered  by  areas  DCIS. 
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Fibroadenoma 

In  Figure  5,  the  one  pre-  and  four  post-contrast  images  of  a  fibroadenoma  are  shown 
as  well  as  the  centroids  and  pixels  associated  with  each  cluster.  Unlike  the  invasive 
carcinoma  in  Fig.  3  and  4,  the  shape  of  the  enhancement  curve  varies  greatly  across  the 
lesion.  This  variation  is  not  apparent  in  the  post-contrast  images  themselves.  However, 
like  the  invasive  carcinoma,  there  is  one  cluster,  3,  which  is  predominantly  the  rim  of  the 
lesion.  On  close  examination  of  the  histology  slides,  this  fibroadenoma  was  found  to  be 
bordered  nearly  completely  by  blood  vessels.  The  high  spatial  variation  of  the 
enhancement  curves  within  the  lesion  may  also  be  a  result  of  the  blood  vessels  inside  the 
lesion. 


Lobular  Carcinoma 

Figure  6  shows  the  centroids  and  membership  map  of  cluster  four  in  12  slices  of  a 
breast  with  extensive  lobular  carcinoma  in  situ  (LCIS).  In  the  post-contrast  images,  this 
breast  would  be  specified  as  having  regional  enhancement.  The  images  how  that  there 
are  two  dominate  regions  of  enhancement,  one  in  the  more  medial  slices  in  the  upper 
region  of  the  breast  and  the  second  in  the  lateral  slices  in  the  middle  of  the  breast.  There 
are  also  blood  vessels  in  the  lateral  slices  at  the  top  of  the  breast  which  proceed  down  the 
surface  of  the  breast  turning  into  the  breast  halfway  down.  This  is  more  easily  seen  in  a 
movie  of  3D  projections.  In  cases  of  regional  enhancement,  it  is  important  to  establish 
connectivity  of  the  patches  of  enhancement  in  order  to  determine  if  it  arises  from  one 
area  and  spreads  of  if  there  are  two  separate  patches  that  developed  independently. 
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Identifying  blood  vessels  is  also  important.  Biopsies  were  obtained  from  the  two  patches 
of  enhancement  and  the  pathologist  identified  the  existence  of  LCIS  in  both  areas. 
Unfortunately,  histology  cannot  determine  whether  these  two  regions  arise  from  the  same 
source  or  are  two  independent  sources.  Based  on  the  3D  FCA  results,  it  appears  as  if 
these  regions  are  connected  spatially. 


Figure  6 


Conclusions 

FCA  was  capable  of  determining  the  location  of  lesions  or  areas  of  LCIS/DCIS. 
However,  by  itself  it  cannot  diagnose  breast  cancer.  Other  parameters  such  as 
architectural  features  and  analysis  of  the  time  course  curves  can  further  aid  in 
distinguishing  benign  from  malignant  disease.  However,  segmentation  using  FCA  can 
point  out  areas  of  concern  to  the  clinician  which  would  aid  in  the  diagnostic  process.  By 
performing  this  analysis  on  a  large  number  of  cases,  uptake  curves  characteristic  of 
various  lesion  types  can  be  identified.  Definition  of  these  curves  is  independent  of  a 
model  for  contrast  uptake.  In  the  future,  the  characteristic  curves  could  be  used  for  both 
diagnosis  and  staging  of  breast  cancer  even  when  the  lesions  are  small. 
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KEY  RESEARCH  ACCOMPLISHMENTS 


•  Determined  optimal  fuzzy  index  and  membership  threshold. 

•  Completed  43  Breast  MRI  studies  which  includes:  benign  lesions, 
ductal  carcinomas,  fibroadenomas,  and  lobular  carcinomas. 
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REPORTABLE  OUTCOMES 


Since  this  is  the  first  year  of  the  project,  publication  of  the  results  detailed  in  this  progress 
have  yet  to  be  published;  two  manuscripts  are  in  progress.  During  this  first  year  period, 
the  graduate  student  (PI)  of  this  research  fellowship  has  published  the  results  of  his 
previous  academic  rotation.  Since  stipend  support  was  received  during  the  period  of  these 
publications,  they  are  appended. 

•  Sodium  and  Proton  MR  Properties  of  Cartilage  During  Compression.  Regatte 
RR.  Kaufman  JH.  Noyszewski  EA.  Reddy  R.  Journal  of  Magnetic  Resonance 
Imaging.  10:961-967,  1999. 

•  A  novel  approach  to  observing  articular  cartilage  deformation  in-vitro  via 
magnetic  resonance  imaging.  Kaufman  JH.  Regatte  RR.  Bolinger  L.  Kneeland 
JB.  Reddy  R.  Leigh  JS.  Journal  of  Magnetic  Resonance  Imaging.  9(5):653- 
662,1999. 

•  Sodium  NMR  evaluation  of  articular  cartilage  degradation.  Insko  EK. 
Kaufman  JH.  Leigh  JS.  Reddy  R.  Magnetic  Resonance  in  Medicine.  41(1):30- 
34,  1999. 
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Original  Research 


A  Novel  Approach  to  Observing  Articular 
Cartilage  Deformation  In  Vitro  via  Magnetic 
Resonance  Imaging 

Jonathan  H.  Kaufman,  MS,*  Ravinder  Reddy  Regatte,  PhD,  Lizann  Bolinger,  PhD, 
J.  Bruce  Kneeland,  MD,  Ravinder  Reddy,  PhD,  and  John  S.  Leigh,  PhD 


The  design  of  a  pressure  cell  that  compresses  a  cartilage 
specimen  in  one  dimension  within  an  imaging  magnet  is 
presented.  One-dimensional  projection  images  in  a  direc¬ 
tion  perpendicular  to  the  articular  surface  of  the  cartilage 
specimen  were  used  to  generate  a  uniaxial  confined  defor¬ 
mation  creep  curve  for  normal  and  trypsin- degraded  carti¬ 
lage  specimens  during  a  continuous  0.690  MPa  (100  psi) 
pressure  application.  The  resulting  curves  are  shown  to  fit 
a  two  time  constant  viscoelastic  model  well  and  also 
indicate  that  the  elastic  modulus  of  cartilage  decreases 
and  the  deformation  rate  increases  upon  trypsin  proteoly¬ 
sis.  Furthermore,  cartilage  permeability  is  shown  as  a 
function  of  cartilage  strain  for  both  the  normal  and  trypsin- 
degraded  case.  Several  two-dimensional  slice-selective  im¬ 
ages  were  collected  both  before  and  after  80  minutes  of 
continuous  compression.  These  images  were  used  to  evalu¬ 
ate  the  relative  changes  in  the  spin-lattice,  Tl,  and  spin- 
spin,  T2,  relaxation  time  constant  maps  for  both  normal 
and  degraded  cartilage  specimens  in  response  to  compres¬ 
sion.  The  results  of  this  study  demonstrate  the  utility  of  a 
novel,  non-magnetic,  cartilage  compression  device  and 
also  support  the  validity  of  a  simple  two-component  rheo¬ 
logical  model  of  articular  cartilage.  J.  Magn.  Reson. 
Imaging  1999;9:653-662.  ©  1999  Wiley-Liss,  Inc. 

Index  terms:  osteoarthritis;  cartilage;  compression;  permeabil¬ 
ity;  viscoelastic;  MRI 

ARTICULAR  CARTILAGE  is  found  at  the  surface  of 
contacting  bones  in  the  synovial  joints  of  mammals.  Its 
role  is  primarily  to  allow  low  friction  joint  mobility.  A 
second  function  of  articular  cartilage  is  to  absorb  me¬ 
chanical  shock  when  synovial  joints  experience  impact 
during  activities  such  as  jogging.  Osteoarthritis  is  a 
condition  in  which  articular  cartilage  degradation  oc¬ 
curs,  leading  to  further  joint  abnormalities  and  in  many 
cases  to  pain  and  loss  of  function  (1).  A  non-invasive 
modality  for  detecting  the  earliest  stages  of  osteoarthri¬ 
tis  has  yet  to  be  developed  (2).  Such  a  non-invasive 
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technique  would  be  clinically  useful  in  assessing  joint 
abnormalities  and  would  also  be  useful  in  assessing  the 
efficacy  of  novel  therapeutics. 

In  its  normal  condition,  articular  cartilage  is  a  remark¬ 
able  tissue  in  the  sense  that  it  is  extremely  firm  while 
being  composed  primarily  of  water  (3).  Approximately 
80%  of  the  mass  of  articular  cartilage  is  water.  The 
remainder  of  the  cartilage  mass  is  composed  for  the 
most  part  of  the  solid  matrix  of  collagen  type  II  mol¬ 
ecules  and  large  proteoglycan  aggregates  (4).  The  com¬ 
position  of  cartilage  is  actively  maintained  by  chondro¬ 
cytes,  although  compared  with  most  tissues,  cartilage  is 
relatively  acellular.  The  firm  nature  of  articular  cartilage 
results  from  a  high  osmotic  pressure.  The  abundant 
proteoglycan  aggregates  are  highly  sulfated,  leading  to  a 
high  negative  fixed  charge  density  within  the  cartilage 
matrix.  This  fixed  charge  density  attracts  a  large  num¬ 
ber  of  cations  (primarily  Na+)  into  the  cartilage  matrix. 
The  resulting  large  number  of  ions  within  the  cartilage 
matrix  attracts  water  at  a  high  pressure  to  satisfy  the 
Donnan  equilibrium  (5) .  However,  there  is  evidence  that 
other  factors  may  contribute  to  the  compressive  strength 
of  cartilage  (6).  This  high  pressure  of  water  within  the 
cartilage  matrix  gives  the  cartilage  its  stiffness,  while 
the  integrity  of  cartilage  is  maintained  by  the  high 
tensile  strength  of  the  collagen  framework. 

The  earliest  stages  of  osteoarthritis  are  often  charac¬ 
terized  by  a  reduction  of  proteoglycan  content  (7-9). 
This  compositional  change  is  thought  to  be  associated 
with  initial  cartilage  softening,  which  is  consistent  with 
the  intuitive  notion  that  a  reduction  in  proteoglycan 
content  would  reduce  the  osmotic  pressure  within  the 
cartilage.  Subsequent  to  this  initial  stage,  the  cartilage 
softening,  in  combination  with  mechanical  wear,  leads 
to  a  disruption  of  the  collagen  framework,  leading  to 
erosion  of  the  articular  surface  of  the  cartilage  followed 
by  fissure  and  cracking. 

The  normal  function  of  articular  cartilage  involves 
compression.  As  articular  cartilage  is  compressed,  the 
fluid  contained  within  the  cartilage  matrix  is  extruded 
into  the  joint  space.  The  extruded  fluid  acts  as  a 
lubricant  between  the  contact  points  of  opposing  bone 
interfaces  (10).  It  is  possible  that  cartilage  that  is  too 
soft  will  compress  to  an  unnatural  extent  and  provide  a 
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mechanical  advantage  for  natural  shear  forces  to  dis¬ 
rupt  the  collagen  framework.  Additionally,  soft  cartilage 
may  not  be  able  to  recover  its  original  shape  quickly 
enough  to  continue  providing  fluid  to  a  bone  interface 
during  repeated  loading.  The  experiment  described  in 
this  manuscript  characterizes  several  magnetic  reso¬ 
nance  imaging  properties  of  articular  cartilage  speci¬ 
mens  while  simultaneously  acquiring  information  re¬ 
garding  deformation. 

MATERIALS  AND  METHODS 
Pressure  Cell  Design 

The  pressure  cell,  shown  in  Fig.  1,  consists  of  a  pneu¬ 
matically  activated  ceramic  piston  confined  in  a  cylindri¬ 
cal  chamber.  The  cartilage  was  positioned  in  the  pres¬ 
sure  cell  with  its  articular  surface  in  contact  with  a 
porous  glass  filter  (Corning  31001-10F).  The  porous 
glass  filter  led  to  a  drain  that  allowed  extruded  cartilage 
fluid  to  leave  the  cell.  The  pressure  cell  design 


inhibits  lateral  cartilage  strain  and  only  permits  fluid  to 
leave  the  cartilage  via  the  articular  surface. 

The  cylinder  and  cap  of  the  pressure  cell  were  made  of 
Delrin  (Dupont).  This  polymer  exhibits  no  proton  NMR 
signal,  is  nonmetallic,  and  has  an  electrical  resistivity 
on  the  order  of  10 14  (1-c m  and  a  dissipation  factor  of 
0.0048  at  106  Hz.  The  piston  was  made  of  a  machinable 
ceramic,  Macor  (Corning).  This  ceramic  has  an  electri¬ 
cal  resistivity  on  the  order  of  1016  fl-cm  and  a  loss 
tangent  of  10-3  at  1  kHz  and  10~4  at  8.5  GHz.  The  pore 
size  of  the  glass  filter  was  chosen  such  that  it  was 
impermeable  to  the  solid  matrix  of  the  cartilage  but 
freely  permeable  to  the  cartilage  fluid.  The  compressed 
air  that  pneumatically  activated  the  piston  was  pre¬ 
vented  from  traveling  around  the  sides  of  the  piston  by  a 
Teflon  seal.  Pressure  was  applied  via  a  compressed  air 
tank  and  was  measured  by  a  digital  pressure  gauge 
(Omega  DPG-500)  on  the  regulator.  A  solenoid  radiofre¬ 
quency  (RF)  coil  was  wrapped  around  the  cylinder  of  the 
pressure  cell.  The  pressure  cell  was  placed  inside  an  1 1 
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Figure  1 .  Cross-sectional  diagram  of  MR-compatible  pressure  cell.  Compressed  air  exerts  a  pressure  on  a  piston  that  is  confined 
in  a  cylindrical  shaft.  The  piston,  when  responding  to  the  pressure  from  the  compressed  air,  pushes  on  the  specimen  that  is 
positioned  between  the  cylinder  and  a  Pyrex  filter,  which  allows  cartilage  fluid  to  extrude  from  the  cell. 
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cm  diameter  local  gradient  set  capable  of  generating  8 
G/cm.  The  cell  and  gradient  set  were  then  placed  in  the 
1  meter  bore  of  a  2  T  Oxford  whole-body  imaging  magnet 
interfaced  to  a  custom-built  NMR  spectrometer. 

Uniaxial  Confined  Deformation  Creep  Experiment 

The  uniaxial  confined  deformation  creep  experiment 
was  chosen  as  a  method  to  measure  the  viscoelastic 
properties  of  the  cartilage  specimens.  This  experiment 
was  first  employed  on  cartilage  specimens  by  Mow  and 
coworkers  (11).  In  this  experiment,  pressure  is  applied 
to  the  cartilage  in  the  form  of  a  step  function.  Prior  to 
pressure  application,  a  series  of  two-dimensional  90  °- 
180  °  spin-echo  images  was  acquired.  Upon  compres¬ 
sion,  one-dimensional  90  -180  °spin-echo  image  projec¬ 
tions  were  acquired  every  5  seconds  for  a  period  of  80 
minutes.  The  compression  was  driven  by  a  pressure  of 
0.690  MPa  (100  psi),  which  was  established  within  the 
first  20  seconds  of  the  one-dimensional  projection- 
image  acquisition  period  and  was  maintained  for  the 
duration  of  the  creep  experiment.  After  the  80  minutes 
of  constant  pressure  application,  the  pressure  of  0.690 
MPa  was  maintained  for  a  subsequent  set  of  two- 
dimensional  image  acquisitions. 

Image  Acquisition 

Slice-selective  images  of  each  cartilage  specimen  in  a 
plane  coincident  with  a  diameter  of  the  cartilage  disc 
and  perpendicular  to  the  cartilage  surface  were  ac¬ 
quired.  Slice  thickness  was  2  mm.  In  each  case  the 
cartilage  specimen  was  oriented  such  that  its  articular 
surface  was  parallel  to  B0,  the  main  magnetic  field.  All 
two-dimensional  images  were  acquired  with  a  standard 
90  -180  °spin-echo  RF  pulse  sequence.  A  10  kHz  band¬ 
width  was  used.  The  orientations  of  frequency  and 
phase  encoding  are  indicated  in  Fig.  1.  Each  image 
employed  a256  X  256  array  of  pixels.  All  images  were 
acquired  using  a  2  T  full-body  imaging  magnet  inte¬ 
grated  with  a  custom-built  spectrometer.  Local  gradi¬ 
ents  sets  that  enclose  the  pressure  cell  were  con¬ 
structed  that  provide  up  to  8  gauss/cm.  All  RF  pulses 
were  applied  via  a  solenoidal  coil  tuned  to  the  proton 
frequency  having  six  turns  and  a  diameter  of  2.5  cm. 
Imaging  data  were  analyzed  using  the  Interactive  Dis¬ 
play  Language  interpreter  on  a  Silicon  Graphics  work¬ 
station. 

Relaxation  Time  Mapping 

We  constructed  maps  of  the  spin-lattice,  Tl,  and  spin- 
spin,  T2,  relaxation  times  for  each  position  of  the 
cartilage  imaging  slice  corresponding  to  an  image  pixel. 
In  general,  the  relaxation  time  maps  were  constructed 
from  fitting  Eqs.  [1]  and  [2] 

I  =  A(1  -  exp  (-  TR/T1))  +  B  (1) 

/  =  C  exp  (-  TE/T2)  (2) 

where  I  is  image  intensity,  TR  is  the  repetition  time,  TE 
is  the  echo  time,  and  A,  B,  C  are  arbitrary  constants,  to  a 


series  of  90  -180  °spin-echo  images  employing  different 
echo  times  and  different  repetition  times.  Five  images 
were  used  to  construct  each  Tl  map.  All  were  acquired 
with  an  echo  time  of  20  msec  and  had  the  following 
repetition  times:  0.5,  1 .0,  1.5,  2.0,  and  2.5  seconds.  Six 
images  were  used  to  construct  each  T2  map.  All  were 
acquired  with  a  repetition  time  of  2.0  seconds  and  had 
the  following  echo  times:  20,  30,  40,  50,  60,  and  70 
msec.  The  collected  images  were  fit  to  Eqs.  [1]  and  [2] 
using  a  gradient  expansion  algorithm  for  a  non-linear 
least  squares  fit.  It  is  important  to  note  that  this 
technique  for  mapping  the  T2  relaxation  time  constant 
is  not  independent  of  diffusion  effects,  and  thus  the  T2 
values  that  we  report  deviate  from  true  T2  values  since 
our  method  has  not  incorporated  a  diffusion-weighting 
factor  (12,13).  As  early  as  1954,  Carr  and  Purcell  (14) 
demonstrated  that  the  contribution  of  diffusion  when 
observing  T2  decay  with  a  series  of  individual  spin  echo 
events  was 

Myd  =  -  M0  exp  [y2Dt3/12n2],  (3) 

where  Myd  is  the  contribution  of  diffusion,  y  is  the 
gyromagnetic  ratio,  D  is  the  diffusion  constant,  t  is 
duration  that  a  spin  ensemble  experiences  a  gradient, 
and  n  is  the  number  of  echoes.  Considering  our  experi¬ 
mental  parameters  of  gradients  no  larger  than  8  G/cm, 
and  using  a  diffusion  constant  of  10-5  cm2/sec,  we 
estimate  that  diffusion  can  influence  our  T2  measure¬ 
ments  by  a  factor  of  approximately  10%. 

In  addition,  given  the  fixed  set  of  repetition  times,  our 
method  for  mapping  Tl  has  an  inherent  limitation  on 
the  range  of  Tl  values.  This  limitation  is  especially 
important  when  considering  Tl  values  near  the  lowest 
repetition  time  of  0.5  seconds.  A  calibration  of  our  Tl 
mapping  technique  with  copper  sulfate-doped  agarose 
phantoms  of  known  T 1  and  cartilage  specimens  (where 
we  measured  bulk  Tl  with  an  inversion  recovery  tech¬ 
nique)  indicated  that  our  Tl  mapping  procedure  yielded 
Tl  values  accurate  to  within  an  approximate  50% 
range. 

Cartilage  Preparation 

All  cartilage  specimens  were  obtained  from  the  patella  of 
bovine  calves  (ages  ranged  from  6  months  to  a  year) 
within  10  hours  from  sacrifice.  We  obtained  a  total  of  10 
specimens,  5  of  which  were  treated  with  trypsin  as 
described  in  the  next  section.  Cartilage  specimens  were 
harvested  from  the  flat  portion  of  the  articular  surface  of 
the  patella  with  a  diamond  core  drill.  With  a  small 
spatula,  each  disc  was  peeled  from  the  cartilage-bone 
interface.  Prior  to  insertion  into  the  pressure  cell,  each 
specimen  was  equilibrated  in  10  mL  of  phosphate- 
buffered  saline  (Sigma  P-4417)  at  room  temperature  for 
30  minutes. 


Trypsin  Degradation 

In  an  effort  to  simulate  articular  cartilage  in  an  early 
arthritic  state,  we  exposed  five  cartilage  specimens  to  a 
moderate  proteolytic  regime.  Since  it  has  been  hypoth- 
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esized  that  early  stages  of  osteoarthritis  involve  the  loss 
of  proteoglycan  content  selectively,  we  chose  a  pro¬ 
teolytic  agent  that  would  selectively  deplete  the  pro¬ 
teoglycan  content  of  the  cartilage  while  leaving  the 
collagen  network  relatively  intact.  The  proteolytic  en¬ 
zyme  trypsin  is  known  to  cleave  the  core  proteins  of  the 
proteoglycans  while  having  little  effect  on  the  type  II 
collagen  fibers  that  comprise  the  collagenous  network  of 
articular  cartilage.  The  digested  cartilage  specimens 
referred  to  in  this  manuscript  were  generated  by  incubat¬ 
ing  each  cartilage  specimen  in  a  10  mL  volume  of  1 
mg/mL  trypsin  (Sigma  T -8253)/ phosphate-buffered  sa¬ 
line  (Sigma  P-44 1 7)  solution  for  7 5  minutes  at  37  G  with 
slight  agitation. 

RESULTS 

Consistency  of  Results  Among  Samples 

The  five  undigested  specimens  had,  with  respect  to  each 
other,  a  variation  of  maximal  strains  of  10%  and  a 
variation  of  deformation  creep  time  constant  of  40%. 
The  results  from  the  five  trypsin-treated  specimens  had 
a  similar  degree  of  variation. 

Viscoelastic  Properties  of  Normal  and  Degraded 
Articular  Cartilage 

One  method  of  measuring  the  mechanical  properties  of 
a  viscoelastic  material  is  with  the  uniaxial  confined 
deformation  creep  experiment.  In  this  method,  a  con¬ 
stant  pressure  is  applied  to  the  viscoelastic  material  as 
the  deformation  of  the  material  is  measured  as  a  func¬ 
tion  of  time  (15).  In  the  context  of  the  experiment  with 
cartilage,  a  first  step  is  to  keep  the  experimental  param¬ 
eters  one-dimensional  such  that  the  analytical  aspects 
are  manageable  (16).  This  is  primarily  possible  with 
cartilage  as  a  consequence  of  the  tissue’s  natural  defor¬ 
mation  mechanism.  Cartilage  can  be  thought  of  as  a 
substance  that,  in  bulk,  is  incompressible  unless  fluid 
extrusion  is  permitted.  We  thus  confined  the  cartilage 
discs  laterally,  while  allowing  fluid  to  extrude  from  the 
articular  surface.  In  this  way,  both  the  pressure  and  the 
deformation  are  contained  in  one  dimension,  the  direc¬ 
tion  perpendicular  to  the  articular  surface. 

To  measure  the  compressive  response  of  the  cartilage 
disc  to  the  applied  pressure,  we  derived  the  thickness  of 
the  cartilage  specimen  from  the  one-dimensional  90  °- 
180  °  spin-echo  image  projections  that  were  acquired 
every  5  seconds  over  a  period  of  80  minutes.  Six  of  these 
projections  are  shown  in  Fig.  2  as  examples.  The  pres¬ 
ence  of  signal  intensity  in  these  projections  indicates 
the  presence  of  protons,  and  for  the  most  part  also 
indicates  the  presence  of  the  cartilage.  Toward  the 
left-hand  side  of  the  peak  is  a  low-intensity  plateau  that 
corresponds  to  the  water  within  the  pores  of  the  perme¬ 
able  glass  filter.  The  representation  of  this  projection 
does  not  have  sharp  edges.  Two  experimental  factors 
contribute  to  this  lack  of  sharpness.  First,  the  cartilage 
specimen  was  not  exactly  perpendicular  to  the  fre¬ 
quency-encoding  direction.  The  tilt  was  slightly  less 
than  3  °.  Second,  a  susceptibility  artifact  existed  at  the 
interface  between  the  cartilage  specimen  and  the  ce¬ 
ramic  piston  as  well  as  between  the  cartilage  specimen 
and  the  glass  filter. 
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Figure  2.  One-dimensional  90  -180  °spin-echo  projection  im¬ 
ages  of  the  cartilage  specimen.  Six  projections  are  shown  at 
various  times  after  onset  of  pressure  application,  as  indicated. 
The  plateau  on  the  left  corresponds  to  water  held  within  the 
glass  filter  shown  in  Fig.  1 .  The  left-hand  slope  of  the  projection 
corresponds  to  the  articular  surface  of  the  cartilage  specimen. 
The  right-hand  side  of  the  projection  corresponds  to  the 
cartilage-bone  boundary.  The  compressive  force  is  from  right 
to  left  against  the  stationary  glass  filter. 

The  thickness  of  the  cartilage  was  estimated  as  the 
full  width  of  the  projection  at  half  its  maximal  value.  The 
thickness  calculation  entailed  a  linear  interpolation  of 
signal  intensity  for  regions  between  those  represented 
by  pixels.  The  thickness  calculation  was  performed 
independently  for  each  one-dimensional  projection.  The 
one-dimensional  projection  images  were  not  collected 
beyond  the  80-minute  data  acquisition  period;  it  is 
likely,  based  on  the  trend  of  the  creep  curves,  that  the 
cartilage  specimens  would  have  compressed  further 
given  more  time. 

We  arrived  at  use  of  the  full  width  at  half- maximum 
for  the  cartilage  thickness  by  starting  with  the  assump¬ 
tion  that  our  signal  acquisition  is  independent  of  posi¬ 
tion  within  the  field  of  view.  We  therefore  convoluted  a 
hypothetical  rectangular  object  in  frequency  space  with 
various  passive  filter  functions  putatively  representing 
the  effect  of  our  measurement  system  (including  rota¬ 
tions  of  the  rectangular  object  that  correspond  to  the 
cartilage  tilt)  and  found  that  the  signal  intensity  at  the 
positions  corresponded  to  the  edges  of  the  initial  box 
function  always  contained  the  value  of  half  the  maxi¬ 
mum  height  of  the  resulting  curve. 

Given  that  we  have  deduced  the  cartilage  thickness  at 
5 -second  intervals  during  the  deformation  creep  experi¬ 
ment,  we  calculated  strain  at  each  of  these  time  points, 
defined  as  the  fractional  change  in  thickness  of  the 
cartilage  from  its  initial  non-compressed  state  as  de¬ 
scribed  by  Eq.  [4] 


where  e,-  is  strain  at  time  i,  Lj  is  cartilage  thickness  at 
time  i,  and  L0  is  the  initial  cartilage  thickness. 

In  doing  so,  we  obtained  the  deformation  creep  curves 
shown  in  Fig.  3,  which  are  similar  to  those  found  in 
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Figure  3.  Deformation  creep  curve  of  normal  and  degraded 
cartilage  specimens.  Strain,  defined  as  the  fractional  change  in 
cartilage  thickness  resulting  from  a  constant  pressure  applica¬ 
tion  of  0.690  MPa,  is  shown  as  a  function  of  time.  The  solid 
lines  represent  biexponential  least  squares  fit  to  the  cartilage 
thickness  calculations,  represented  as  circles.  The  curve  for 
the  normal  cartilage  specimen  is  labeled  n,  and  the  curve 
corresponding  to  the  trypsin-degraded  cartilage  specimen  is 
labeled  d.  The  standard  deviations  of  strain  measurements  are 
shown  as  dotted  lines. 

previous  purely  mechanical  studies  (17,18).  Most  nota¬ 
bly,  Mow  et  al  (18)  observed  a  15%  strain  for  an  applied 
pressure  of  0.1  MPa.  The  shape  of  the  creep  curve 
generated  in  Mow  et  al’s  experiment  is  similar  to  the 
shape  we  obtained  and  has  a  similar  time  constant.  We 
have  fit  these  curves  to  a  general  exponential  decay 
function  given  by  Eq.  [5], 

e  =  njj  |l  -  exp  j—  “|  +  |l  —  exp  (-  TJ  +  (5) 

where  e  is  strain,  tis  time,  and  (irq,  m2,  m3,  m4,  m5)  are  fit 
parameters.  Two  time  constants  were  required  for  an 
adequate  fit.  The  least  squares  fit  results  are  given  in 
Table  1 . 

Effects  of  Compression  on  Imaging  Parameters 

Spin-lattice,  Tl,  and  spin-spin,  T2  magnetic  resonance 
relaxation  time  constant  maps  were  generated  from  a 
series  of  slice-selective  spin-echo  images  acquired  both 
before  and  after  compression  of  the  cartilage.  As  dis- 


Table  1 

Deformation  Creep  Curves  Fit  to  a  Biexponential  Decay* 


Normal  specimen 

Degraded  specimen 

m\ 

0.18397 

0.27486 

rrh 

43.256 

29.206 

rrh 

0.030851 

0.080153 

m4 

0.19517 

0.16205 

ms 

3.7087 

1.5768 

*The  deformation  creep  data  of  Fig.  3  were  fit  to  the  biexponential 
decay  function  given  in  Eq.  [5]. 


cussed  in  the  methods  section,  these  measurements  of 
Tl  and  T2  are  used  for  comparative  purposes.  The 
values  reported  may  contain  inaccuracies  since  the 
measurement  methods  were  non-standard  and  were 
not  calibrated.  These  maps,  along  with  proton  density- 
weighted  images,  are  shown  in  Fig.  4  for  representative 
normal  and  trypsin-degraded  cartilage  specimens.  In 
the  cases  of  both  normal  and  degraded  cartilage  speci¬ 
mens,  the  proton  density-weighted  images  displayed  a 
fairly  uniform  image  intensity  prior  to  and  subsequent 
to  compression.  In  contrast,  the  relaxation  time  maps 
had  a  uniform  intensity  only  after  compression.  Prior  to 
compression,  Tl  was  fairly  heterogeneous  throughout 
the  imaging  slice,  and  high  values  of  T2  were  predomi¬ 
nantly  located  in  the  region  of  the  articular  surface  in 
both  the  normal  and  trypsin- degraded  cases. 

Since  the  relaxation  time  constants  throughout  the 
imaging  slice  of  the  cartilage  specimens  were  not  uni¬ 
form  prior  to  compression,  we  have  calculated  Tl  and 
T2  histograms  of  the  normal  and  degraded  cartilage 
specimens  prior  to  and  after  compression.  Tl  histo¬ 
grams  for  both  specimens  before  and  after  compression 
are  shown  in  Fig.  5.  We  should  note  that  experimental 
factors  could  have  contributed  to  the  relaxation  time 
distributions.  For  example,  the  width  of  the  Tl  histo¬ 
gram  could  represent  noise  in  the  Tl  map  rather  than 
an  actual  Tl  distribution,  especially  since  Tl  in  carti¬ 
lage  has  been  shown  to  be  relatively  homogeneous  (19). 
When  testing  this  notion  by  evaluating  the  pixel-by¬ 
pixel  reproducibility  of  Tl  maps  of  cartilage  specimens 
and  agarose  phantoms,  we  found  relatively  little  consis¬ 
tency.  However,  the  fact  that  the  Tl  distribution  nar¬ 
rows  upon  compression  could  implicate  the  existence  of 
Tl  heterogeneity.  Although  there  was  a  reduction  in  Tl 
for  both  specimens  upon  compression,  the  change  in  Tl 
for  the  degraded  specimen,  upon  compression,  was 
substantially  more  dramatic  than  that  corresponding  to 
the  normal  specimen. 

T2  histograms  for  both  specimens  before  and  after 
compression  are  shown  in  Fig.  6.  Unlike  the  case  for  Tl , 
the  reductions  in  T2  for  both  the  normal  and  degraded 
specimens  were  similar.  However,  T2  measurements  for 
the  degraded  specimen  both  before  and  after  compres¬ 
sion  were  slightly  higher  than  those  for  the  normal 
specimen.  As  is  indicated  in  the  images  of  Fig.  4,  the 
histograms  of  Figs.  5  and  6  likewise  indicate  that  both 
relaxation  time  constants  become  more  uniform  upon 
compression  for  normal  and  degraded  cartilage. 

From  the  Tl  and  T2  histograms,  we  have  calculated 
average  values  of  Tl  and  T2  prior  to  and  after  compres¬ 
sion  for  both  the  normal  and  degraded  specimens. 
These  values  are  given  in  Table  2. 

Permeability  Calculation 

The  data  from  the  deformation  creep  experiment  per¬ 
mits  us  to  calculate  the  permeability  of  the  cartilage 
assuming  that  the  deformations  are  within  the  elastic 
limit.  The  most  direct  approach  to  calculating  permeabil¬ 
ity  within  this  experimental  context  is  to  use  Darcy’s 
law,  which  states  that  the  velocity  of  the  fluid  in  a 
permeable  solid  matrix  is  linearly  proportional  to  the 
pressure  gradient  (17).  The  permeability  of  the  matrix, 
k,  is  defined  as  the  constant  of  proportionality.  This 
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Figure  4.  Images  and  relaxation  time  maps  of  normal  and  trypsin-degraded  cartilage  specimens.  Each  slice  is  along  a  diameter  of 
the  cartilage  specimen.  The  left  side  of  each  image  corresponds  the  region  of  cartilage  adjacent  to  the  ceramic  piston  of  the 
pressure  cell.  The  right  side  of  each  image  corresponds  to  the  articular  surface  of  the  cartilage  adjacent  to  the  glass  filter.  Image 
intensities  of  the  normal  and  degraded  spin-echo  images  are  shown  with  the  same  image  intensity  scaling  before  and  after 
compression,  respectively.  All  T1  maps  have  the  same  scaling,  and  likewise,  all  T2  maps  have  the  same  scaling. 


relationship  is  expressed  symbolically  in  Eq.  [6]. 

v  =  kVp  (6) 

where  v  is  the  velocity,  and  p  is  the  pressure.  All  three 
variables  in  this  relationship  are  assumed  to  be  func¬ 
tions  of  both  space  and  time. 

Since  our  experimental  setup  is  designed  to  study  the 
deformation  of  cartilage  in  one  dimension  (ie,  the  carti¬ 
lage  is  confined  such  that  it  cannot  expand  laterally) ,  we 
can  treat  Eq.  [6]  in  a  one-dimensional  fashion  where  x  is 
position.  Since  our  experiment  collects  data  at  discrete 
time  points,  it  is  useful  to  include  time  explicitly  in  a 
discrete  fashion,  as  indicated  by  Eq.  [7]  where  the 
subscript  i  indicates  a  particular  point  in  time: 


14 


It  will  be  useful  to  express  Eq.  [7]  completely  in  terms  of 
parameters  derived  directly  from  a  deformation  creep 
experiment,  namely,  L0,  the  original  thickness  of  the 
cartilage;  P,  the  constant  pressure  applied  to  the  carti¬ 
lage;  Lf,  the  equilibrium  thickness  of  the  cartilage  under 
pressure  P;  At,  the  time  step  between  each  data  acquisi¬ 
tion;  and  ei#  the  strain  at  time  step  i. 


Since  the  fluid  can  be  considered  incompressible 
for  our  purposes,  for  the  sake  of  simplicity,  we  make 
the  approximation  that  the  velocity  of  the  fluid  at  all 
points  within  the  cartilage  is  uniform  and  can  be  ex¬ 
pressed  as 


v,  = 


AL, 

At 


(8) 


where 


AL/  =  U  ~  L/-!  (9) 

and  Li  is  the  thickness  of  the  cartilage  at  time  i.  Further¬ 
more,  since  our  working  definition  of  strain  is 


e,  = 


L/  L0 

L0 


(10) 


a  substitution  of  Eq.  [10]  into  Eq.  [8]  yields  an  explicit 
expression  for  y- 


V/  = 


L0Ae/ 


At 
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Figure  5.  Apparent  spin-lattice  relaxation  time  constant,  Tl,  histograms  of  the  normal  and  degraded  cartilage  specimens  both 
before  and  after  compression. 


where 


equilibrium  strain  for  a  given  pressure 


Ae,  =  e,  -  e,.,.  (12) 

Let  us  assume  for  simplicity  that  the  elastic  properties 
of  the  cartilage  specimen  are  homogeneous.  With  this 
assumption,  and  since  there  are  no  pressure  sources  or 
sinks  within  the  cartilage  volume,  the  pressure  gradient 
throughout  the  cartilage  volume  is  constant.  Thus,  the 
pressure  can  be  expressed  as  a  linear  function  of  space 

p  =  Ax  4-  B  (13) 


E  = 


P 

{Lf  -  L0)/L0  ’ 


(16) 


Substituting  Eq.  [16]  into  Eq.  [15]  and  dividing  by  the 
thickness  of  the  cartilage  at  time  i  leads  to  the  uniform 
pressure  gradient,  given  by 


/dp 

[dx 


=  (P/Q  i  - 


ei^0  \ 

(Lf~L0)  ]’ 


(17) 


where  A  and  B  are  defined  by  the  boundary  conditions. 
Not  all  of  this  pressure,  however,  affects  the  velocity  of 
the  interstitial  fluid,  since  a  portion,  pe,  is  absorbed  by 
the  elastic  resistance  of  the  cartilage 


Substituting  Eqs.  [17]  and  [11]  into  Eq.  [7],  we  obtain  an 
explicit  expression  for  permeability  in  experimentally 
derived  terms 


pe  =  Ee  (14) 

where  E  is  the  elastic  modulus  of  the  cartilage.  The 
change  in  pressure  across  the  thickness  of  the  cartilage 
is  thus 


A  p  =  P-pe=P-Ee  (15) 

where  P  is  the  constant  applied  pressure.  The  elastic 
modulus,  E,  of  the  cartilage  can  be  defined  by  the 


k:  = 


LoAe,- 


/PAA 

o 

<lT 

U) 

'“"a 

1 

i 

— < 

(18) 


The  biexponential  fits  to  the  deformation  creep  data, 
when  substituted  into  Eq.  [18],  provide  the  plots  of 
permeability  as  a  function  of  strain  given  in  Fig.  7.  The 
figure  indicates  that,  in  the  uncompressed  state,  the 
trypsin-degraded  specimen  is  approximately  twice  as 
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Figure  6.  Apparent  spin-spin  relaxation  time  constant,  T2,  histograms  of  the  normal  and  degraded  cartilage  specimens  both 
before  and  after  compression. 


permeable  to  the  normal  specimen  and  that  the  in¬ 
crease  in  permeability  with  respect  to  strain  is  sharper 
than  that  of  the  normal  specimen.  The  permeability  of 
the  specimen  in  the  region  of  high  strain  (greater  than 
0.3)  is  unlikely  to  be  accurately  represented  by  this 
method  since  the  fluid  velocity  in  this  range  becomes 
small.  Additionally,  in  the  region  of  high  strain,  our 
experimental  conditions  exceed  the  elastic  limit  for 
articular  cartilage  as  reported  by  Mow  et  al  (18).  Figure 
8  shows  the  permeability,  calculated  in  the  same  man- 


Table  2 

Changes  in  Apparent  Relaxation  Parameters  With  Compression* 


T1  (sec) 

T2  (msec) 

Normal  before 

0.99  ±  0.38 

66.4  ±  47.2 

Normal  after 

0.62  ±  0.55 

38.4  ±  35.6 

Degraded  before 

0.91  ±  0.36 

71 .9  ±46.5 

Degraded  after 

0.60  ±  0.57 

46.9  ±  39.7 

*The  averages  of  spin-latice  (T1)  and  spin-spin  (T2)  relaxation  time 
constants  were  calculated  from  the  maps  shown  in  Fig.  4. 


ner,  as  a  function  of  time  for  the  deformation  of  both  the 
normal  and  degraded  specimen.  This  figure  shows  that 
during  the  course  of  deformation,  the  degraded  speci¬ 
men  is  initially  more  permeable  than  the  normal  speci¬ 
men,  but  after  approximately  5  minutes  of  sustained 
compression,  the  degraded  specimen  loses  its  permeabil¬ 
ity  with  respect  to  the  normal  specimen. 

DISCUSSION 
MR  Properties 

The  experiments  described  in  this  manuscript  expand 
on  earlier  progress  to  identify  early  manifestations  of 
osteoarthritis  using  MR  imaging  such  as  observations  of 
MR  signal  intensity  changes  across  the  depth  of  the 
cartilage  (20)  and  measurements  of  MR  relaxation  time 
constants  of  normal  and  degraded  cartilage  (21).  MR 
relaxation  times  have  also  previously  been  observed  to 
change  in  a  depth-dependent  fashion  in  normal  articu¬ 
lar  cartilage  specimens  during  compression  (22).  Our 
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Figure  7.  Permeability  of  the  normal  and  degraded  cartilage  specimens  as  a  function  of  strain,  defined  as  the  absolute  value  of 
the  fractional  change  in  cartilage  thickness  during  compression. 


present  observations  of  the  MR  properties  of  bovine 
specimens  of  articular  cartilage  are  consistent  with 
these  previous  studies.  However,  we  were  unable  to 
observe  consistent  depth-dependent  results  with  all 
specimens.  This  inconsistency  could  result  from  the 
likelihood  that  juvenile  cartilage  specimens  (which  we 
used)  are  likely  to  be  structurally  distinct  from  adult 
cartilage  (3).  In  addition,  with  our  technique,  it  was  difficult 
to  assess  the  MR  characteristics  of  the  cartilage  at  the 


Figure  8.  Permeability  of  the  normal  and  degraded  cartilage 
specimens  as  a  function  of  time. 


articular  surface  since  some  imaging  artifacts  did  arise 
from  the  interface  of  the  cartilage  with  the  glass  filter. 

Viscoelastic  Properties 

There  are  a  multitude  of  theories  describing  the  me¬ 
chanical  properties  of  articular  cartilage,  including  the 
electrostatic  interactions  among  negatively  charged  pro¬ 
teoglycans  (23),  the  viscous  interaction  between  the 
cartilage  fluid  and  its  solid  matrix  (11),  and  the  micro- 
structural  organization  of  cartilage  (24).  It  has  been 
proposed  that  purely  viscoelastic  models  of  articular 
cartilage  may  not  be  adequate  for  a  comprehensive 
treatment  of  the  mechanical  properties  of  cartilage  (25). 
However,  the  distinction  between  poroelastic  and  vis¬ 
coelastic  deformations  could  be  a  matter  of  semantics. 
Viscoelastic  models  are  found  to  be  consistent  with 
numerous  experimental  observations  (26). 

We  chose  to  apply  a  simple  rheological  model  to  our 
data  and  expand  it  as  necessary.  Since  the  data  from 
our  deformation  creep  experiment  at  first  glance  approxi¬ 
mate  a  function  that  exponentially  rises  to  a  fixed  value, 
we  chose  the  Voigt  (15)  viscoelastic  model  shown  in  Fig. 
9,  which  is  described  by  Eq.  [19]. 


where  e  is  strain,  P  is  pressure,  E  is  the  elastic  modulus, 
v  is  the  viscosity,  and  c  is  an  arbitrary  constant. 
However,  we  found  that  two  time  constants  were  neces¬ 
sary  for  an  adequate  fit  to  our  observations.  Thus,  the 
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Figure  9.  Voigt  viscoelastic  model.  This  basic  viscoelastic 
model  consists  of  an  elastic  spring  arranged  in  parallel  to  a 
viscous  dashpot.  As  indicated,  E  is  the  elastic  modulus,  and  v 
is  the  viscosity. 


mathematical  description  of  two  Voigt  solids  in  series 
adequately  fit  the  data  as  shown  in  Fig.  3. 

Our  experiment  was  designed  to  mimic  physiological 
pressures  that  are  well  beyond  the  elastic  limit  of 
cartilage,  which  does  place  limits  on  the  present  physi¬ 
cal  interpretation  of  a  viscoelastic  model.  However,  we 
can  state  that  the  deformation  creep  observations  are 
consistent  with  the  notion  that  the  cartilage  rheology 
can  be  described  by  a  linked  deformation  of  two  compo¬ 
nents. 

In  general  the  observed  changes  in  permeability  and 
elastic  moduli  resulting  from  degradation  are  consis¬ 
tent  with  mechanical  studies  involving  alternative  ar¬ 
thritic  models,  such  as  a  canine  surgically  induced 
arthritic  model  (27),  as  well  as  other  enzymatic  degrada¬ 
tion  models  (28).  Furthermore,  our  observations  that 
permeability  increases  with  trypsin  digestion  is  consis¬ 
tent  with  the  notion  that  a  reduction  in  fixed  charge 
density  within  the  cartilage  matrix  induces  an  increase 
in  permeability  (29).  We  have  described  an  approach  to 
investigating  cartilage  functionality  in  vitro  that  could 
lend  insight  into  the  design  of  experiments  in  vivo. 
Furthermore,  the  degree  of  consistency  between  the 
experimental  results  and  a  two  time  constant  viscoelas¬ 
tic  model  lends  credence  to  the  notion  that  the  histologi¬ 
cal  zones  of  articular  cartilage  have  distinct  mechanical 
properties. 
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Sodium  and  Proton  MR  Properties  of  Cartilage 
During  Compression 

Ravinder  Reddy  Regatte,  PhD,*  Jonathan  H.  Kaufman,  MS,  Elizabeth  A.  Noyszewski,  PhD, 
and  Ravinder  Reddy,  PhD 


Proton  and  sodium  MR  relaxation  times  of  bovine  articular 
cartilage  specimens  were  measured  as  a  function  of  proteo¬ 
glycan  (PG)  depletion  and  as  a  function  of  mechanical 
compression.  Proton  and  sodium  relaxation  times  of  nor¬ 
mal  cartilage  were  compared  with  relaxation  times  of 
PG-depleted  cartilage  to  evaluate  the  significance  of  PG 
depletion-induced  changes  in  MR  relaxation  parameters. 
These  comparisons  were  conducted  for  both  uncompressed 
and  mechanically  compressed  states.  The  mechanical  com¬ 
pressions  were  performed  with  an  MR-compatible  pressure 
cell  and  evaluated  dynamically  via  interleaved  one-dimen¬ 
sional  proton  and  sodium  MR  projection  imaging.  The 
comparisons  indicate  that  sodium  relaxation  parameters 
are  sensitive  to  PG  depletion  when  cartilage  is  in  a  mechani¬ 
cally  compressed  state  or  an  uncompressed  state.  In  con¬ 
trast,  proton  relaxation  parameters  do  not  change  signifi¬ 
cantly  with  PG  depletion  when  cartilage  is  in  an 
uncompressed  state.  However,  during  mechanical  compres¬ 
sion,  proton  T2  becomes  sensitive  to  PG  depletion.  These 
results  support  the  potential  of  sodium  magnetic  reso¬ 
nance  imaging  (MRI)  as  a  possible  modality  for  obtaining 
imaging  contrast  related  to  PG  depletion.  The  results  also 
indicate  the  potential  of  proton  MRI  to  provide  such  con¬ 
trast  if  the  image  acquisition  is  conducted  in  conjunction 
with  a  mechanical  compression  via  physical  exercise.  J. 
Magn.  Reson  Imaging  10:961-967,  1999. 

©  1999  Wiley-Liss,  Inc. 
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ARTICULAR  CARTILAGE  is  a  dense  connective  tissue 
composed  primarily  of  water,  collagen  fiber  bundles, 
and  proteoglycan  aggregates.  Cartilage  functions  as  the 
load-bearing  surface  of  all  synovial  joints.  The  normal 
function  of  articular  cartilage  involves  compression.  As 
cartilage  is  compressed,  the  fluid  contained  within  the 
matrix  is  extruded  into  the  joint  space.  This  extruded 
fluid  acts  as  a  lubricant  between  the  contact  points  of 
opposing  bone  interfaces  (1).  The  disease  process  of 
osteoarthritis  (OA)  involves  the  degeneration  of  articular 
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cartilage.  The  earliest  putative  changes  in  cartilage  due 
to  OA  are  a)  softening  of  cartilage,  b)  a  decrease  in 
proteoglycan  content,  c)  an  increase  in  water  content,  d) 
changes  in  cartilage  thickness,  and  e)  an  increase  in 
permeability  (2-5).  Each  of  these  changes  has  been 
associated  with  an  alteration  in  the  mechanical  proper¬ 
ties  of  cartilage  (6-8).  Therefore,  there  has  been  substan¬ 
tial  interest  in  studying  the  mechanical  properties  of 
normal  and  degenerated  articular  cartilage. 

Several  investigators  have  used  proton  MRI  to  observe 
cartilage  image  dynamics  during  compression  (9,10). 
However,  they  could  only  report  qualitative  effects  be¬ 
cause  MR  relaxation  time  constants  were  not  mea¬ 
sured.  Previously,  we  have  estimated  the  relative 
changes  in  proton  MR  relaxation  parameters  of  articu¬ 
lar  cartilage  specimens  as  a  function  of  compression 
(11).  This  previous  work  indicated  that  both  proton  T1 
and  T2  decrease  with  compression.  However,  this  indi¬ 
cation  was  based  on  non-standard  methods  for  measur¬ 
ing  relaxation  parameters.  In  the  present  work  we  use 
standard  inversion  recovery  and  Hahn  spin-echo  tech¬ 
niques  to  obtain  robust  values  for  proton  relaxation  param¬ 
eters  of  cartilage  during  mechanical  compression. 

In  addition  to  proton,  sodium  is  an  important  nucleus 
in  cartilage  physiology.  Several  investigators  have  mea¬ 
sured  the  MR  properties  of  sodium  as  a  function  of 
cartilage  degeneration  (3,12-15).  Such  studies  with 
sodium  have  yet  to  incorporate  the  effects  of  mechanical 
compression.  Therefore,  in  addition  to  the  proton  mea¬ 
surements,  we  have  extended  the  objective  of  the  pre¬ 
sent  study  to  obtain  measurements  of  sodium  relax¬ 
ation  parameters  as  a  function  of  compression.  In  both 
cases,  we  examined  the  effects  of  enzymatic  proteogly¬ 
can  depletion  on  the  observed  compression-induced 
trends  in  relaxation  times. 


MATERIALS  AND  METHODS 
Cartilage  Specimens 

Ten  1  -cm-diameter  cylindrical  cartilage  plugs  were  har¬ 
vested,  each  from  the  flat  surface  of  a  bovine  patella.  All 
patellae  were  stored  at  20 °C  prior  to  the  experiment.  It  is 
unknown  what  the  effect  of  freezing  is  on  cartilage 
specimens.  Each  cartilage  specimen  was  sectioned  with 
a  diamond  core  drill  and  then  peeled  from  the  cartilage- 
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bone  interface  with  a  spatula.  The  specimens  were 
subsequently  inspected  to  ensure  that  there  was  no 
visible  damage  and  that  all  cartilage  layers  remained 
intact.  Each  specimen  was  equilibrated  in  phosphate- 
buffered  saline,  pH  7.4  (Sigma,  St.  Louis,  MO)  for  30 
minutes. 

Five  of  the  cartilage  specimens  were  exposed  to  enzy¬ 
matic  digestion.  Since  it  has  been  hypothesized  that 
early  stages  of  OA  involve  the  selective  loss  of  proteogly¬ 
can  (PG)  content,  we  chose  a  proteolytic  agent  that 
would  deplete  the  proteoglycan  content  of  the  cartilage 
while  leaving  the  collagen  network  relatively  intact.  The 
digested  cartilage  specimens  were  generated  by  incubat¬ 
ing  each  cartilage  specimen  in  5  ml  of  1 .0  mg  trypsin/ ml 
phosphate-buffered  saline  solution  (Sigma)  for  60  min¬ 
utes.  The  end-stage  digestion  medium  was  analyzed  for 
proteoglycan  content  via  a  spectrophotometric  assay 
(16).  In  each  case,  the  PG  concentration  in  the  media 
was  approximately  1 .3  mg/ml,  which  corresponds  to  an 
approximate  30%  PG  depletion. 

Pressure  Cell  Design  and  Construction 

The  pressure  cell  consists  of  a  pneumatically  activated 
ceramic  piston  and  porous  glass  filter  in  a  cylindrical 
Delrin  chamber,  as  shown  in  Fig.  1 .  The  details  of  the 
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Figure  1.  Cross-sectional  diagram  of  MR-compatible  pressure 
cell.  Compressed  air  exerts  a  pressure  on  a  piston  that  is 
confined  in  a  cylindrical  shaft.  The  piston,  when  responding  to 
the  pressure  from  the  compressed  air,  pushes  on  the  specimen 
that  is  positioned  between  the  cylinder  and  a  Pyrex  filter,  which 
allows  cartilage  fluid  to  extrude  from  the  cell. 


pressure  cell  are  given  elsewhere  (11).  Each  cartilage 
plug  was  placed  between  the  piston  and  filter  such  that 
the  articular  surface  of  the  cartilage  was  adjacent  to  the 
porous  glass  filter.  The  pressure  cell  was  oriented  such 
that  the  cartilage  articular  surface  was  parallel  to  the 
main  magnetic  field  (B0).  The  ceramic  piston  com¬ 
pressed  each  cartilage  specimen,  while  the  porous  glass 
filter  allowed  the  extrusion  of  cartilage  fluid.  A  double- 
tuned  solenoid  radiofrequency  (RF)  coil  was  incorpo¬ 
rated  into  the  pressure  cell.  The  transmit/receive  sole- 
noidal  RF  coil  was  double-tuned  to  the  Larmor 
frequencies  of  H-l  (86.1  MHz)  and  Na-23  (22.7  MHz)  as 
described  by  Schnall  et  al  (17).  The  pressure  cell  was 
placed  inside  an  1 1  -cm-diameter  local  gradient  set 
capable  of  generating  8  G/cm.  The  cell  and  gradient  set 
were  then  placed  in  the  center  of  a  2  T -Oxford  whole- 
body  imaging  magnet  interfaced  to  custom-built  NMR 
spectrometer.  The  pressure  cell  was  activated  pneumati¬ 
cally,  via  a  regulated  (0.69  MPa)  compressed  air  source. 
The  pressure  was  measured  at  the  source  with  a  digital 
pressure  gauge  (Omega  DPG-500). 


Proton  and  Sodium  Relaxation  Time 
Meas  uremen  ts 

The  standard  inversion  recovery  technique  was  em¬ 
ployed  to  obtain  proton  T1  measurements.  For  normal 
cartilage  specimens,  the  inversion  time  was  varied  by 
step  times  of  200  msec  before  compression  and  100 
msec  during  compression.  For  the  PG-depleted  carti¬ 
lage  specimens,  the  inversion  time  varied  by  step  times 
of  250  msec  before  compression  and  125  msec  during 
compression.  In  each  case  30  free  induction  decays 
(FIDs)  were  acquired.  Each  FID  contained  1024  data 
points  at  a  sampling  rate  of  0. 1  msec. 

To  measure  proton  T2,  the  standard  Hahn  spin-echo 
technique  (18)  was  employed.  For  normal  cartilage,  the 
echo  time  (tau)  was  varied  by  25  msec  before  compres¬ 
sion  and  10  msec  during  compression.  For  PG-depleted 
cartilage,  tau  was  varied  by  20  msec  before  compression 
and  8  msec  during  compression.  Spectra  at  30  time 
points  were  acquired  for  each  T2  measurement.  For 
each  time  point,  200  FIDs  were  averaged  to  obtain  a 
spectrum. 

Although  sodium  exhibits  bi- exponential  T 1  in  biologi¬ 
cal  tissues,  due  to  the  small  difference  between  the  two 
T Is,  it  is  difficult  to  measure  each  independently.  There¬ 
fore,  measurements  of  sodium  spin  lattice  relaxation 
time  constants,  Tl,  were  carried  out  using  the  same 
standard  inversion-recovery  methods  used  for  proton. 
However,  since  the  Tl  of  sodium  is  several  orders  of 
magnitude  shorter  than  that  of  proton,  the  inversion 
time  was  varied  by  a  step  time  of  6  msec  before  compres¬ 
sion  and  3  msec  during  compression  for  both  normal 
and  PG-depleted  cartilage  specimens.  Thirty  inversion 
time  points  were  acquired  for  each  Tl  measurement; 
200  FIDs  were  acquired  and  averaged  for  each  inversion 
time.  Each  FID  contained  1024  data  points  at  a  sam¬ 
pling  rate  of  0. 1  msec. 

The  transverse  relaxation  of  sodium  is  bi-exponential 
and  thus  has  two  time  constants:  T2f  and  T2s.  It  has 
been  shown,  using  triple  quantum  filtered  spectra,  that 
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only  the  slow  component,  T2s,  is  sensitive  to  compres¬ 
sion  (19).  Therefore  we  used  the  same  standard  Hahn 
spin-echo  technique  for  the  measurement  of  sodium 
transverse  relaxation  times  as  was  used  as  in  the  case 
for  proton  transverse  relaxation  times.  For  normal  carti¬ 
lage,  tau  was  varied  by  4  msec  before  compression  and  2 
msec  during  compression.  For  PG-depleted  cartilage, 
tau  was  varied  by  3  msec  before  compression  and  1 
msec  during  compression.  Spectra  at  30  time  points 
were  acquired  for  each  T2  measurement.  For  each  time 
point,  200  FIDs  were  averaged  to  obtain  a  spectrum. 

Projection  Images 

One-dimensional  proton  and  sodium  projection  images, 
in  a  direction  perpendicular  to  the  articular  surface, 
were  acquired  with  the  pulsed  gradient  spin-echo  pulse 
sequence  (PGSE) .  The  proton  imaging  parameters  were 
as  follows:  TR  1  second,  TE  20  msec,  number  of  excita¬ 
tions  (NEX)  1.  For  sodium  projections,  the  imaging 
parameters  were  as  follows:  TR  100  msec,  TE  3  msec, 
NEX  800.  The  acquisitions  of  proton  and  sodium  projec¬ 
tion  images  were  interleaved.  The  total  acquisition  time 
for  one  set  of  interleaved  proton  and  sodium  projection 
images  was  81  seconds  (1  second  for  the  proton  projec¬ 
tion  image  and  80  seconds  for  the  sodium  projection 
image).  One  hundred  interleaved  projection  images  were 
collected  for  a  total  duration  of  approximately  135 
minutes. 

RESULTS 

Demonstration  of  Compression 

Figure  2  shows  a  series  of  one-dimensional  proton 
projection  images  of  cartilage  acquired  before  and  dur¬ 
ing  the  deformation.  These  projections  were  acquired 
every  81  seconds  for  a  period  of  135  minutes.  The 
ceramic  piston  compressed  the  cartilage  from  the  right 
side  toward  the  filter  on  the  left.  The  low-intensity 
plateau  on  the  left  side  of  the  projection  image  corre¬ 
sponds  to  the  position  of  the  glass  filter.  The  projection 
images  narrow  with  the  pressure  onset,  and  the  signal 
intensity  drops  as  a  function  of  time. 

Figure  3  shows  a  series  of  one-dimensional  sodium 
projection  images  of  cartilage  before  and  during  defor¬ 
mation.  As  in  the  case  of  proton  projections,  the  sodium 
projections  were  acquired  every  8 1  seconds  for  a  period 
of  135  minutes.  With  pressure  application,  the  sodium 
projection  images  narrow,  and  sodium  signal  de¬ 
creases. 


Figure  2.  One-dimensional  90-180°  spin-echo  proton  projec¬ 
tion  images  of  (a)  a  normal  cartilage  specimen,  and  (b)  a 
PG-depleted  cartilage  specimen.  Five  projections  are  shown  at 
various  times  after  onset  of  pressure  application.  The  compres¬ 
sive  force  is  from  right  to  left  against  the  stationary  glass  filter. 
The  left-hand  slope  of  the  projection  corresponds  to  the  articu¬ 
lar  surface  of  the  cartilage  specimen.  The  right-hand  slope  of 
the  projection  corresponds  to  the  cartilage-bone  boundary. 

regions  between  those  represented  by  pixels.  Strain  can 
be  defined  in  the  following  three  ways: 


Evidence  of  Enzymatic  Induction  of  Cartilage 
Softening 

Figure  4  shows  the  cartilage  strain  obtained  from  the 
proton  and  sodium  projection  images  as  a  function  of 
compression  duration  for  both  normal  and  PG-depleted 
cartilage.  The  strain  calculations  employed  the  consid¬ 
eration  that  the  thickness  of  a  cartilage  specimen  is 
approximately  the  full  width  at  half  maximum  of  the 
corresponding  one-dimensional  projection  images 
shown  in  Figs.  2  and  3.  The  thickness  calculation 
employed  linear  interpolation  of  signal  intensity  for 


(Lp-LJ 

0/2)(Lo  +  LJ  ’ 


(2) 

(3) 


where  L0  is  the  original  thickness,  and  Lx  is  the  thick¬ 
ness  after  compression.  In  this  manuscript,  we  use  the 
definition  of  strain  expressed  by  Eq.  [1].  Cartilage  strain 
was  evaluated  at  8 1  second  intervals. 
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(b)  PG  depleted 


Figure  3.  One-dimensional  spin-echo  sodium  projection  im¬ 
ages  of  (a)  a  normal  cartilage  specimen,  and  (b)  a  PG-depleted 
cartilage  specimen.  Five  projections  are  shown  at  indicated 
times  after  the  onset  of  pressure  application.  The  specimen 
orientation  is  identical  to  that  of  Fig.  2. 

At  the  onset  of  pressure  application,  cartilage  strain 
rapidly  increased  to  a  steady-state  condition.  The  PG- 
depleted  cartilage  specimens  were  systematically  more 
compliant  than  the  normal  specimens.  The  experimen¬ 
tal  creep  data  for  both  sodium  and  proton  were  fitted  to 
a  simple  bi-exponential  decay  with  two  time  constants 

given  by 

e  =  A„  +  Aj|l  -  exp  j-Tjj  +  A2|l  -  exp  j-TjJ  (4) 

where  e  is  strain,  t  is  time,  and  Aq,  Alt  A2,  0i  and  02  are  fit 
parameters  (given  in  Table  1).  This  model  mimics  the 
mechanical  processes  of  articular  cartilage  during  con¬ 
fined  compression  as  shown  previously  (11). 

The  shapes  of  the  sodium  and  proton  creep  curves  are 
similar.  The  PG-depleted  specimens  were  almost  twice 
as  mechanically  compliant  as  the  normal  specimens. 


The  lower  signal-to-noise  ratios  of  the  sodium  projec¬ 
tions  compared  with  the  proton  projections  is  reflected 
as  the  large  standard  error  of  the  sodium  creep  curves. 

Effects  of  Compression  on  MR  Relaxation  Times 

Figure  5  shows  typical  curve  fits  when  calculating 
relaxation  time  constants.  Table  2  presents  a  summary 
of  the  data  obtained  regarding  the  effects  of  compres¬ 
sion  and  PG  depletion  on  cartilage  proton  and  sodium 
MR  relaxation  times.  The  four  inherent  parameters: 
proton-Tl,  proton-T2,  sodium-Tl,  and  sodium-T2  are 
delineated  for  the  four  physio-mechanical  cartilage 
states:  normal-uncompressed,  normal-compressed,  PG- 
depleted-uncompressed,  and  PG-depleted-compressed. 

Each  of  the  four  parameters  were  evaluated,  using  a 
two-tailed  t- test,  to  determine  whether  PG  depletion  had 
a  significant  effect  on  any  of  the  MR  time  constants.  The 
resulting  P  values  are  given  in  the  last  column  of  Table 
2.  In  the  uncompressed  state,  the  sodium  MR  time 
constants  show  a  much  stronger  PG  depletion-induced 
distinction  than  do  the  proton  MR  time  constants. 
There  is  no  evidence  that  proton-Tl  changes  with  PG 
depletion,  and  scarce  evidence,  if  any,  that  proton-T2 
changes.  On  the  other  hand,  both  sodium-Tl  and 
sodium-T2  increase  with  PG  depletion  to  a  statistically 
significant  extent. 

Compressing  cartilage  provides  a  means  of  introduc¬ 
ing  a  PG  depletion-induced  change  in  proton-T2.  In  the 
compressed  state  (using  the  same  pressure),  the  pro¬ 
ton-!^  of  the  PG-depleted  cartilage  is  shorter  than  that 
of  normal  cartilage.  Furthermore,  the  change  in  pro- 
ton-T2  associated  with  compression  is  substantially 
more  pronounced  for  PG-depleted  cartilage  than  it  is  for 
normal  cartilage. 

The  contribution  of  compression  in  distinguishing  the 
effects  of  PG  depletion  on  sodium  MR  relaxation  param¬ 
eters  is  less  pronounced;  the  change  in  sodium-T2, 
associated  with  compression,  is  similar  for  the  normal 
and  PG-depleted  specimens. 

The  data  also  provide  an  evaluation  of  the  changes  in 
proton  and  sodium  relaxation  time  constants  associ¬ 
ated  with  compression.  To  compliment  the  values  dis¬ 
played  in  Table  2,  we  report  the  statistical  significance 
of  the  compression-associated  changes  in  relaxation 
time  in  Table  3. 

DISCUSSION 

The  central  aim  of  this  study  was  to  identify  PG  deple¬ 
tion-induced  trends  in  proton  and  sodium  MR  relax¬ 
ation  times.  These  trends  could  serve  as  potential 
contrast  mechanisms  that  indicate  early  cartilage  degen¬ 
eration  and  thus  early  signs  of  osteoarthritis.  Since  the 
normal  functionality  of  articular  cartilage  involves  com¬ 
pression,  we  have  measured  the  MR  relaxation  times  of 
articular  cartilage  in  both  uncompressed  and  mechani¬ 
cally  compressed  states.  One  motivation  for  the  intro¬ 
duction  of  mechanical  compression  is  that  it  has  been 


Figure  4.  Typical  relaxation  time  constant  data  fit  to  exponen¬ 
tial  functions:  (a)  proton T1 ,  (b)  sodium T1 ,  (c,  d)  proton  T2,  and 
(e,  f)  sodium  T2. 
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Table  1 

Deformation  Creep  Curve  Fit  Parameters  of  Fig.  4  Using  Eq.  [4] 


Nucleus 

PG 

depleted 

Ao 

A, 

a2 

9i 

02 

Proton 

No 

0.0364 

0.175 

0.159 

60.025 

8.99 

Proton 

Yes 

0.0385 

0.127 

0.127 

18.153 

28.09 

Sodium 

No 

0.0046 

0.157 

0.155 

8.182 

54.042 

Sodium 

Yes 

-0.1171 

0.395 

0.395 

3.66 

3.67 

hypothesized  that  cartilage  softening  may  be  the  earli¬ 
est  known  indicator  of  osteoarthritis. 

The  method  of  this  study  employed  one-dimensional 
projection  imaging  to  monitor  and  verify  cartilage  com¬ 
pression.  Furthermore,  the  degree  of  compression  and 
the  rate  at  which  compression  occurred  for  a  given 
compressive  load  were  calculated  via  the  construction 
of  deformation  creep  curves.  These  curves  permitted 


(a)  From  proton  data 


(b)From  sodium  data 


Figure  5.  Deformation  creep  curve  of  a  normal  and  a  PG- 
depleted  cartilage  specimen  obtained  from  proton  and  sodium 
projections.  The  curve  for  the  normal  cartilage  specimen  is 
labeled  (n),  and  the  curve  corresponding  to  the  PG-depleted 
cartilage  specimen  is  labeled  (d).  The  solid  line  represent 
bi-exponential  least  squares  fits  of  the  raw  data  to  Eq.  [41,  and 
the  dotted  lines  correspond  to  standard  deviations  of  strain 
measurements.  - 


Table  2 

Sodium  and  Proton  Relaxation  Parameters* 


Full 

compression 

Normal 

SD 

PG 

depleted 

SD 

P 

value 

Proton  T1  (sec) 

No 

1.54  0.06 

1.52 

0.09  0.658 

Yes 

0.99  0.10 

1.07 

0.09  0.217 

(%  change) 

-35.6 

7.9 

-29.5 

5.8 

0.205 

Proton  T2  (msec) 

No 

82.4 

2.7 

90.6 

10.6 

0.155 

Yes 

38.8 

2.6 

29.6 

6.2 

0.028 

(%  change) 

-52.9 

2.5 

-67.1 

6.8 

0.007 

Sodium  T1  (msec) 

No 

24.0 

1.6 

30.6 

2.7 

0.003 

Yes 

14.8 

1.9 

24.8 

3.4 

0.001 

(%  change) 

-38.2 

8.2 

-19.1 

7.0 

0.004 

Sodium  T2  (msec) 

No 

18.0 

3.2 

25.7 

3.8 

0.008 

Yes 

11.2 

0.8 

15.5 

2.1 

0.009 

(%  change) 

-36.6 

8.6 

-39.3 

4.1 

0.548 

*Relaxation  times  are  given  for  proton  T 1 ,  proton  T2,  sodium  T 1 ,  and 
sodium  T2  for  cartilage  in  the  normal  and  mechanically  compressed 
state.  A  comparison  is  drawn  between  these  parameters  for  normal 
cartilage  and  cartilage  that  has  been  enzymatically  treated  to  induce 
proteoglycan  depletion.  Standard  deviations  of  all  relaxation  times 
are  given  in  the  adjacent  columns  labeled  SD.  In  each  case,  five 
cartilage  specimens  were  used.  The  relaxation  times  reported  are  the 
corresponding  means.  In  addition  to  relaxation  times,  the  percent 
change  between  the  normal  and  compressed  state  is  given.  For  each 
normal  vs.  PG-depleted  comparison,  the  result  of  a  two-tailed  t- test  is 
given  in  the  rightmost  column,  labeled  P,  the  value  in  this  column 
represents  the  probability  that  the  normal  and  PG-depleted  values 
could  have  arisen  from  sampling  a  single  population. 

the  verification  that  enzymatic  PG  depletion  effectively 
increased  cartilage  softness.  The  establishment  of  carti¬ 
lage  compression  and  PG  depletion  permitted  an  evalu¬ 
ation  of  PG  depletion-induced  changes  in  proton  and 
sodium  MR  characteristics  of  cartilage  in  the  normal 
state  and  during  mechanical  compression. 

The  results  of  the  present  study  have  implications  for 
clinical  research  in  cartilage  imaging.  Standard  clinical 
proton  imaging  that  employs  basic  T 1  and  T2  contrast 
mechanisms  is  unlikely  to  distinguish  PG  depletion  in 
articular  cartilage  since  we  have  demonstrated  that 
proton  T1  and  T2  do  not  vary  significantly  with  PG 
depletion.  However,  we  have  demonstrated  that,  when 
cartilage  is  in  a  mechanically  compressed  state,  a  signifi¬ 
cant  change  in  proton  T2  does  occur  with  PG  depletion. 
This  result  indicates  that  it  may  be  possible  to  identify 
regions  of  articular  cartilage  that  have  a  low  PG  content 
with  T2 -weighted  proton  imaging,  if  the  cartilage  of 
interest  is  in  a  compressed  state.  A  limitation  worth 
mentioning  is  that  compressed-state  proton  T2  changes 


Table  3 

Statistical  Significance  of  the  Change  in  Relaxation  Time  Constants 
Associated  With  Compression 


Nuclei 

Normal/degraded 

Relaxation 
time  constant 

Fvalue 

Proton 

Normal 

T1 

1.60E-05 

Normal 

T2 

5.00E-09 

Degraded 

T1 

4.00E-05 

Degraded 

T2 

3.20E-05 

Sodium 

Normal 

T1 

3.50E-05 

Normal 

T2 

0.0056 

Degraded 

T1 

0.018 

Degraded 

T2 

0.002 
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associated  with  PG  depletion  are  on  the  order  of  magni¬ 
tude  of  their  standard  deviation  among  samples.  There¬ 
fore,  relaxation  time  constant,  alone,  may  provide  infor¬ 
mation  on  relative  PG  depletion  (ie,  a  measure  of  cartilage 
degradation  over  time)  but  will  require  additional  infor¬ 
mation  to  make  an  absolute  prediction  of  PG  depletion. 
In  vivo,  mechanical  compression  of  articular  cartilage 
can  be  induced  physiologically  by  having  subjects  per¬ 
form  certain  exercises  prior  to  imaging.  Experimental 
evidence  supports  the  notion  that  the  rate  of  cartilage 
recovery  is  slow  enough  to  permit  MRI  measurements  of 
cartilage  thickness  dynamics  during  recovery  (20). 

The  results  also  support  the  notion  that  sodium 
imaging  could  have  a  substantial  impact  on  cartilage 
evaluation  in  a  clinical  setting.  We  have  demonstrated 
that  both  sodium  T1  and  T2  change  significantly  as  a 
function  of  PG  depletion  whether  or  not  the  cartilage  is 
in  a  compressed  state.  One  caveat  to  this  result  is  that 
sodium  T1  and  T2  both  increase  with  PG  depletion,  and 
thus  each  may  offset  the  effect  of  the  other  depending  on 
the  specifics  of  an  imaging  pulse  sequence.  The  sodium 
relaxation  parameters  that  we  report  can  aid  in  the 
development  of  an  optimal  sodium  pulse  sequence  for 
imaging  cartilage  PG  loss. 

In  addition  to  clinical  implications,  the  results  re¬ 
ported  in  this  study  are  relevant  to  understanding  the 
physiological  mechanisms  of  cartilage  deformation  and 
NMR  relaxation  mechanisms  in  cartilage.  For  example, 
the  proton  relaxation  time  constants  reported  here  are 
consistent  with  those  of  previous  work  in  which  carti¬ 
lage  compression  was  performed  via  an  osmotic  tech¬ 
nique  prior  to  MR  evaluation  at  atmospheric  pressure 
(21).  The  consistency  of  the  present  results  with  those 
cited  indicates  that  changes  in  compression-induced 
proton  MR  characteristics  can  be  generated  by  both 
mechanical  or  osmotic  pressure  gradients.  This  conclu¬ 
sion  has  implications  for  in  vitro  proton  MR  cartilage 
studies  in  that  investigators  should  be  careful  to  match 
the  osmolarity  of  cartilage  tissue  culture  media  with 
that  of  synovial  fluid.  A  further  implication  is  that  it  is 
possible  that  baseline  measurements  of  cartilage  proton 
relaxation  times  in  vivo  can  be  influenced  by  the  osmo¬ 
larity  by  a  subject’s  synovial  fluid. 

With  the  assumption  that  PG  depletion  is  the  primary 
cause  of  cartilage  softening,  these  results  indicate  that 
it  may  be  possible  to  infer  cartilage  softening  via  an 
indirect  measurement  of  PG  content.  As  mentioned 
above,  there  is  no  direct  way  of  measuring  PG  changes 
using  proton  MR  without  resorting  to  compression. 
However,  using  sodium  MR,  one  can  indirectly  infer 
changes  in  PG  content,  without  resorting  to  compres¬ 
sion.  Two  limitations  of  using  sodium  MRI  are  that  it 
has  a  lower  signal-to-noise  ratio  and  that  it  requires 
stronger  imaging  gradients  to  resolve  thin  objects  such 
as  cartilage.  However,  the  shorter  T1  of  sodium  can  be 
exploited  in  signal  averaging  to  improve  signal-to-noise 
ratio.  High-quality  sodium  images  of  articular  cartilage 
have  been  acquired  in  vivo  in  a  research  setting  (12), 
which  suggests  that  the  technology  can  be  moved  to¬ 
ward  clinical  applications.  This  inference  can  then  be 
related  to  cartilage  softening  and  in  turn  to  a  non- 


compressive  assessment  of  cartilage  functionality.  Fur¬ 
ther  work  in  this  direction  is  in  progress. 
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Sodium  NMR  Evaluation  of  Articular 
Cartilage  Degradation 
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One  of  the  first  effects  of  degenerative  osteoarthritis  is  the  loss 
of  proteoglycans  from  the  matrix  of  articular  cartilage.  Using  a 
model  of  osteoarthritic  change  where  the  cartilage  has  been 
enzymatically  degraded  with  trypsin,  the  sodium  NMR  character¬ 
istics  of  the  cartilage  were  determined  as  a  function  of  changes 
in  the  proteoglycan  content.  The  results  demonstrate  that  the 
single  quantum  sodium  signal  decreases  as  the  proteoglycan 
content  of  the  cartilage  matrix  decreases.  In  addition,  the 
relaxation  characteristics  of  the  sodium  change  such  that  the  7} 
is  longer,  the  7^s  is  longer,  and  the  T2 f  is  shorter.  Short  echo-time, 
Tj -weighted  sodium  images  are  presented  which  demonstrate 
that  this  information  may  be  utilized  to  detect  the  loss  of 
proteoglycans  from  articular  cartilage.  Magn  Reson  Med 
41:30-34, 1999.  ©  1999  Wiley-Liss,  Inc. 

Key  words:  sodium;  cartilage;  osteoarthritis 

INTRODUCTION 

Articular  cartilage  functions  as  the  load-bearing  surface  of 
all  joints.  The  load-bearing  capability  is  a  function  of  the 
proteoglycan  content  of  the  cartilage  matrix.  Negatively 
charged  proteoglycans  in  the  matrix  of  the  cartilage  in¬ 
crease  the  concentration  of  sodium  ions  present  in  articu¬ 
lar  cartilage  to  levels  that  may  approach  300  m M  or  higher 
in  healthy  cartilage  (1,  2).  The  osmotic  pressure  of  the 
sodium  ions  provides  the  cartilage  with  its  load-bearing 
capability. 

In  response  to  IL-1,  chondrocytes  release  metalloprotein- 
ases  that  degrade  proteoglycans  in  the  cartilage  (3-7) .  One 
of  the  earliest  effects  of  degenerative  osteoarthritis  is  the 
loss  of  proteoglycans  from  articular  cartilage  (5-7).  When 
the  proteoglycan  component  of  the  cartilage  decreases, 
there  is  a  simultaneous  decrease  in  the  sodium  concentra¬ 
tion  in  the  cartilage  (1,  8).  In  addition,  the  loss  of  the 
negatively  charged  proteoglycans  from  the  cartilage  may 
alter  the  relaxation  characteristics  of  the  sodium.  Such 
changes  have  been  demonstrated  by  sodium  NMR  (9-1 1). 

The  purpose  of  this  paper  is  to  measure  the  sodium 
spectroscopic  parameters  of  articular  cartilage  as  a  func¬ 
tion  of  changes  in  the  proteoglycan  content  of  the  cartilage 
and  to  determine  which  parameters  change  after  the  loss  of 
even  small  amounts  of  proteoglycans  from  the  matrix  of  the 
cartilage.  The  change  in  the  single  quantum  sodium  signal 
as  a  function  of  changes  in  the  proteoglycan  content  of 
cartilage  is  demonstrated.  The  change  in  the  T\  and  both 
components  of  the  sodium  Tz,  TW  and  T2S,  as  a  function  of 
changes  in  the  proteoglycan  content  of  cartilage  are  also 
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presented.  Sodium  images  are  presented  which  demon¬ 
strate  that  the  loss  of  proteoglycans  from  the  matrix  of  the 
cartilage  may  be  detected  with  single  quantum,  7\- 
weighted,  short  echo-time  methods. 

MATERIALS  AND  METHODS 

Cartilage  Specimens 

Specimens  of  articular  cartilage  were  cored  from  fresh 
bovine  patella  using  a  diamond-tipped  cylindrical  core 
drill  bit  (Accuratus  Ceramics,  Washington,  NJ).  The  result¬ 
ing  disk-shaped  specimen  was  roughly  1  cm  in  diameter. 
Nine  cartilage  specimens  were  obtained  for  evaluation. 
Eight  of  these  specimens  were  dissected  away  from  the 
underlying  bone  and  marrow,  leaving  an  approximately 
2-mm  thick  layer  of  cartilage.  These  eight  specimens 
ranged  in  wet  weight  from  0.23  to  0.33  g  and  were  retained 
for  sodium  spectroscopic  evaluation:  four  control  speci¬ 
mens  and  four  specimens  for  degradation.  One  specimen, 
with  intact  underlying  bone  and  marrow,  was  cut  into  a 
rectangle  with  a  scalpel  blade.  This  specimen  was  2  mm 
wide,  4  mm  thick,  and  1  cm  long.  This  sample  was  then 
divided  in  two  pieces  of  5  mm  length  each.  These  two 
specimens  were  retained  for  sodium  imaging  evaluation. 

Enzymatic  Degradation 

Proteolytic  degradation  of  the  articular  cartilage  was  per¬ 
formed  in  a  trypsin-containing  bath.  Trypsin  degrades 
proteoglycans  by  cleaving  the  core  protein  into  succes¬ 
sively  smaller  units  (12).  The  negatively  charged  gly- 
cosaminoglycans,  side  chains  of  the  core  protein,  are 
liberated  from  the  cartilage  in  this  process.  The  degrada¬ 
tion  bath  was  prepared  with  a  solution  of  0.2  mg/ml  of 
trypsin  (Sigma  Chemical  Co.,  St.  Louis,  MO)  in  a  phos¬ 
phate  buffered  saline  solution  (150  mM  sodium  chloride) 
that  had  been  adjusted  to  pH  7.6.  A  control  solution 
containing  no  trypsin  was  also  prepared.  Five  cartilage 
specimens  were  placed  in  the  trypsin  bath  for  2-h  intervals 
at  30°C  and  0.5  Hz  agitation.  Another  five  specimens  were 
placed  in  the  control  solution  for  2-h  intervals  under 
equivalent  conditions.  After  each  2-h  interval  in  the  bath, 
the  solutions  were  withdrawn  for  spectrophotometric  evalu¬ 
ation  of  glycosaminoglycan  content  (13)  and  the  cartilage 
was  evaluated  with  sodium  NMR:  either  spectroscopy  or 
imaging.  After  data  acquisition,  the  cartilage  was  returned 
to  a  fresh  solution,  either  trypsin  or  control,  for  another  2-h 
interval.  Five  cycles  of  data  acquisition  were  performed  for 
the  specimens  that  were  evaluated  spectroscopically.  These 
specimens  underwent  a  total  of  10  h  of  enzymatic  degrada¬ 
tion  during  the  experiment.  Four  cycles  of  data  acquisition 
were  performed  for  the  specimens  that  were  evaluated  by 
single  quantum  sodium  imaging.  This  specimen  under- 
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went  a  total  of  8  h  of  enzymatic  degradation  during  the 
experiment. 

Spectrophotometric  Evaluation 

Spectrophotometric  evaluation  of  the  glycosaminoglycan 
content  in  the  degradation  solution  was  performed  by 
adding  serial  dilutions  of  the  solution  to  a  5-ml  solution  of 
16  mg/liter  dimethylmethylene  blue  in  a  buffer  of  5 
ml/liter  ethanol,  2  g/liter  sodium  formate,  and  2  ml/liter 
formic  acid  (13).  The  glycosaminoglycan  content  of  the 
degradation  solutions  were  determined  by  observing  the 
change  in  spectrophotometric  absorption  at  535  nm  and 
comparing  the  change  with  that  observed  in  a  set  of 
standard  solutions  with  known  concentrations  of  chondroi- 
tin  sulfate  (Sigma  Chemical  Co.,  St.  Louis,  MO).  The  total 
glycosaminoglycan  content  of  each  piece  of  cartilage  was 
determined  by  enzymatically  degrading  each  piece  for  48  h 
in  multiple  trypsin  solutions,  to  control  for  enzymatic 
autodigestion,  after  all  NMR  data  were  gathered.  After  this 
period  of  time,  no  additional  glycosaminoglycan  could  be 
extracted  from  any  of  the  cartilage  specimens. 

Sodium  Spectroscopy 

All  NMR  experiments  were  performed  in  a  2  T  whole  body 
magnet  (Oxford  Magnet  Technology,  UK)  that  had  been 
interfaced  to  a  custom-built  spectrometer  hardware  and 
software  package.  Each  cartilage  specimen  retained  for 
spectroscopic  evaluation  was  placed  at  the  center  of  a 
custom-built  5-turn,  1-cm  inside  diameter  solenoid  radio¬ 
frequency  coil.  The  coil  was  tuned  to  22.84  MHz,  the 
resonance  frequency  of  sodium-23  at  2  T,  using  a  variable 
capacitor  (Crane,  New  Rochelle,  NY)  and  fixed  capacitors 
(American  Technical  Ceramics,  Huntington  Station,  NY). 

Sodium  spectroscopic  evaluation  of  the  articular  carti¬ 
lage  was  performed  to  evaluate  the  single  quantum  sodium 
signal,  effective  Ti,  and  the  two  components  of  the  biexpo¬ 
nential  TZ\TZ f  and  T2s  via  a  triple  quantum  method.  The 
single  quantum  signal  was  obtained  from  each  specimen 
with  a  single  90°  pulse  experiment  that  averaged  600 
acquisitions  at  a  repetition  time  of  200  ms.  The  triple 
quantum  sodium  signal  was  obtained  from  each  specimen 
with  a  90o-180o-90o-90°  pulse  sequence  that  employed  a 
six-step  phase  cycling  scheme  (14).  The  preparation  time, 
the  time  between  the  first  two  90°  pulses  in  the  sequence, 
was  varied  in  steps  from  0. 1  to  30  ms  to  obtain  antiphase 
biexponential  decay  data  that  were  subsequently  fit  to 
determine  T2f  and  T2s.  Each  triple  quantum  data  point 
averaged  600  total  acquisitions  at  a  repetition  time  of  200 
ms.  The  T{  of  the  sodium  was  determined  with  a  180°-90° 
inversion  recovery  sequence.  The  inversion  time,  the  time 
between  the  initial  180°  and  the  90°  observation  pulse,  was 
varied  in  2-ms  increments  from  an  initial  value  of  2  ms  to  a 
maximal  value  of  100  ms.  Each  data  point  was  signal 
averaged  25  times  at  a  repetition  rate  of  200  ms. 

Sodium  Imaging 

Cartilage  specimens  were  placed  into  a  thin-walled  glass 
tube  (Wilmad  Glass  Co.,  Buena,  NJ)  that  was  positioned  at 
the  center  of  a  custom-built  5-turn,  5-mm  inside  diameter 
solenoid  radiofrequency  coil  positioned  at  the  center  of  a 


custom-built  gradient  set  (15).  The  coil  was  tuned  as 
described  previously. 

Sodium  images  were  obtained  with  a  non-slice  selective 
spin-echo  sequence  described  in  detail  elsewhere  (16).  The 
slice  thickness  was  defined  by  the  thickness  of  the  sample, 
which  was  2  mm.  The  imaging  matrix  acquired  128  phase 
encoding  steps  and  256  data  points  in  each  acquisition. 
The  pulse  sequence  acquired  an  asymmetric  echo  (17)  and 
was  reconstructed  via  a  homodyne  reconstruction  tech¬ 
nique  (18,  19).  The  FOV  was  3  cm  in  the  phase  encoding 
direction  and  1  cm  in  the  frequency  encoding  direction, 
which  was  through  the  thickness  of  the  cartilage.  The  echo 
time  was  1  ms.  The  repetition  time  was  30  ms.  A  total  of  50 
signal  averages  were  acquired  in  3  min  and  12  s. 

Data  Analysis 

The  spectroscopic  sodium  signal  strengths  were  determined 
using  a  principle  component  analysis  algorithm  (20,  21). 
The  Ti  of  the  sodium  was  determined  by  fitting  the  inversion 
recovery  sodium  signal  as  a  function  of  the  time  between 
inversion  and  observation  pulses  using  an  adaptation  of  a 
nonlinear,  least  squares  fit  algorithm  (22)  to  the  function 

S(t)  =  A{  1  -  2e-|/T0  +  B  [1] 

In  this  function,  t  is  the  time  between  the  180°  inversion 
pulse  and  the  observation  pulse  for  the  data  observed  at 
each  time  point,  and  A  and  B  represent  an  amplitude  and 
offset  level,  respectively. 

The  Tz  values  for  the  sodium,  T2f  and  T2s,  were  deter¬ 
mined  by  fitting  the  triple  quantum  sodium  signal  as  a 
function  of  the  preparation  time  in  the  acquisition  se¬ 
quence.  This  yields  a  data  set  that  varies  in  time  as  the 
difference  of  two  decaying  exponentials,  one  of  short-time 
constant,  Tz f,  and  one  of  long-time  constant,  TZs.  The  data 
were  fit,  using  an  adaptation  of  a  nonlinear,  least  squares  fit 
algorithm  (22)  to  the  function 

S(t)  =  A(e~t/T 2s  —  e_i/T2f)  [2] 

where  t  was  the  preparation  time  in  the  triple  quantum 
pulse  sequence  (14)  and  A  is  the  amplitude. 

RESULTS 

The  cartilage  specimens  investigated  in  this  experiment 
were  degraded  from  their  original  proteoglycan  content  down 
to  between  50%  and  70%  of  their  original  content.  The 
single  quantum  signal  and  T2f  both  decreased  with  proteo¬ 
glycan  degradation,  whereas  the  Tx  and  TZs  both  increased 
with  proteoglycan  degradation  of  the  cartilage  specimens. 

The  single  quantum  sodium  signal  as  a  function  of  the 
proteoglycan  content  of  the  cartilage  is  shown  in  Fig.  1. 
The  data  demonstrate  that  the  single  quantum  sodium 
signal  decreased  approximately  in  proportion  to  the  proteo¬ 
glycan  content  of  the  specimens. 

The  two  components  of  the  sodium  T2,  T2f  and  T2s,  were 
determined  from  fitting  the  amplitude  of  the  envelope  of 
the  triple  quantum  sodium  signal  as  a  function  of  the 
preparation  time  in  the  triple  quantum  pulse  sequence. 
Before  any  enzymatic  degradation,  the  T2f  of  the  cartilage 
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SQC  and  Proteoglycan  Content 
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FIG.  1.  The  normalized  SQC  (single  quantum  coherence)  sodium 
signal  in  the  enzymatically  degraded  cartilage  specimens  as  a 
function  of  the  measured  proteoglycan  content  remaining  in  the 
cartilage  after  each  interval  of  enzymatic  degradation. 

specimens  was  1.9  ±  0.1  ms.  The  T2s  before  any  degrada¬ 
tion  was  7.2  ±  0.4  ms.  With  these  values  and  a  delay 
between  the  center  of  the  RF  excitation  pulse  and  the  first 
data  point  of  100  ms,  the  sodium  visibility  was  approxi¬ 
mately  94%.  The  data  in  Fig.  2  demonstrate  that  the  fast 
component  of  the  sodium  T2,  T2 r,  reduced  slowly  with  the 
initial  proteoglycan  degradation  and  then  somewhat  more 
rapidly  as  the  proteoglycan  degradation  increased.  The 
data  in  Fig.  3,  however,  demonstrate  that  the  slower 
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FIG.  2.  The  T2\  of  the  sodium  signal  in  the  enzymatically  degraded 
cartilage  specimens  as  a  function  of  the  measured  proteoglycan 
content  remaining  in  the  cartilage  after  each  interval  of  enzymatic 
degradation. 
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FIG.  3.  The  T2s  of  the  sodium  signal  in  the  enzymatically  degraded 
cartilage  specimens  as  a  function  of  the  measured  proteoglycan 
content  remaining  in  the  cartilage  after  each  interval  of  enzymatic 
degradation. 

component  of  the  sodium  T2,  T2si  was  prolonged  as  the 
proteoglycan  content  of  the  cartilage  decreased. 

Before  any  enzymatic  degradation,  the  T\  of  the  cartilage 
specimens  was  17.3  ±  0.6  ms.  The  Ti  of  sodium  as  a 
function  of  the  proteoglycan  content  of  articular  cartilage  is 
shown  in  Fig.  4.  The  data  show  that  the  Tx  increased  as  the 
proteoglycan  content  of  the  cartilage  decreased.  The  change 
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FIG.  4.  The  7}  of  the  sodium  signal  in  the  enzymatically  degraded 
cartilage  specimens  as  a  function  of  the  measured  proteoglycan 
content  remaining  in  the  cartilage  after  each  interval  of  enzymatic 
degradation. 
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FIG.  5.  Single  quantum  sodium  images  of  a  section  of  articular  cartilage  that  has  been  divided  into  two  pieces.  The  top  piece  in  each  image 
has  been  enzymatically  degraded  in  a  trypsin-containing  bath  at  2-h  intervals.  The  approximate  percentage  of  the  proteoglycan  content 
remaining  in  the  top  piece  is  listed  above  each  image.  The  image  repetition  time  was  30  ms  and  the  echo  time  was  1  ms. 


in  the  Tx  was  clear  even  after  a  small  degree  of  cartilage 
degradation. 

Single  quantum  sodium  images  of  the  two  cartilage 
specimens  retained  for  imaging  are  shown  in  Fig.  5.  The 
image  repetition  time  and  echo  time  were  chosen  to 
maximize  the  decrease  in  signal  of  the  enzymatically 
degraded  specimen  while  still  allowing  for  rapid  signal 
averaging  so  that  high  signal-to-noise  ratio  images  can  be 
produced.  The  images  demonstrate  that  the  sodium  image 
intensity  decreased  significantly  in  the  specimen  that  was 
enzymatically  degraded.  A  region-of-image  analysis  of  200 
pixels  from  the  central  portion  of  the  cartilage  was  per¬ 
formed  to  evaluate  changes  in  image  intensity.  The  results, 
shown  in  Fig.  6,  give  the  image  intensity  as  a  function  of 
changes  in  the  proteoglycan  content  of  the  cartilage. 

The  loss  of  signal  in  the  images  shown  in  Fig.  6  is  a  result 
of  the  change  in  both  the  sodium  concentration  and  the 
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FIG.  6.  The  normalized  signal  intensity  in  a  200-pixel  region  of 
interest  in  the  central  portion  of  the  image  of  the  enzymatically 
degraded  cartilage  specimens  as  a  function  of  the  measured  proteo¬ 
glycan  content  remaining  in  the  cartilage  after  each  interval  of 
enzymatic  degradation.  The  error  bars  are  the  standard  deviation  of 
the  measurements. 


relaxation  characteristics  of  the  sodium.  After  8  h  of 
degradation,  approximately  31%  of  the  proteoglycan  con¬ 
tent  has  been  removed  and  the  sodium  signal  has  de¬ 
creased  by  approximately  45%.  With  a  1-ms  echo  time,  a 
reduction  in  T2f  from  1.9  to  1.4  ms  would  decrease  the 
image  intensity  by  21%,  and  an  increase  in  T2s  from  7.2  to 
10  ms  would  increase  image  intensity  by  4%.  The  exact 
contribution  of  these  two  effects  is  unknown,  however,  as 
the  proportion  of  the  sodium  in  the  bound  and  free  pool  is 
unknown.  An  increase  in  the  T\  from  17.3  to  24  ms,  with  a 
30-ms  repetition  time  and  a  90°  excitation  pulse,  would 
decrease  image  intensity  by  15%.  A  change  in  sodium 
concentration,  which  is  due  to  a  change  in  proteoglycan 
concentration,  would  decrease  the  image  intensity  by 
approximately  31%  at  this  degree  of  degradation.  A  combi¬ 
nation  of  these  effects  is  required  to  explain  the  45% 
decrease  in  image  intensity  observed  in  this  specimen  after 
8  h  of  degradation. 


CONCLUSIONS 

In  this  work,  we  have  demonstrated  the  change  in  sodium 
spectroscopic  parameters  that  occur  as  a  function  of  changes 
in  the  proteoglycan  content  of  articular  cartilage.  The  data 
showed  that  the  single  quantum  sodium  signal  decreased 
approximately  in  proportion  to  the  proteoglycan  content  of 
the  cartilage.  The  single  quantum  sodium  signal  is  known 
to  change  stoichiometrically  with  the  proteoglycan  content 
of  the  cartilage  (1). 

The  change  in  the  sodium  Tx  as  a  function  of  changes  in 
the  proteoglycan  content  was  demonstrated  for  the  first 
time.  These  data  revealed  that  a  prolongation  of  the  Tx 
occurs  with  loss  of  the  proteoglycans  from  the  cartilage  (10, 
1 1)  and  that  the  prolongation  of  the  Tx  occurs  even  with  a 
small  loss  of  the  proteoglycan  content  of  the  cartilage. 
Also,  the  data  revealed  a  change  in  the  transverse  relax¬ 
ation  times,  T2(  (10)  and  T2s,  with  loss  of  proteoglycans 
from  the  cartilage.  The  changes  were  demonstrated  as  a 
function  of  changes  in  the  proteoglycan  content  of  the 
cartilage  for  the  first  time.  These  data  revealed  that  prolon¬ 
gation  in  T2s  also  occurs  relatively  early  in  the  degradation 
process.  The  T2f,  alternately,  became  shorter  with  the  loss 
of  proteoglycans  from  the  matrix  of  the  cartilage.  This 
change,  however,  was  only  clear  with  greater  levels  of 
proteoglycan  loss  from  the  cartilage. 
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These  data  demonstrated  that  certain  characteristics  of 
the  sodium  NMR  signal  may  be  exploited  to  detect  the  loss 
of  the  proteoglycan  component  of  articular  cartilage  matrix 
via  sodium  imaging.  Most  of  the  change  in  the  image  signal 
likely  arises  simply  from  a  decrease  in  the  sodium  content 
of  the  cartilage  secondary  to  a  loss  of  the  proteoglycan 
content  of  the  cartilage.  However,  the  change  in  the  relax¬ 
ation  characteristics  of  sodium  as  a  function  of  the  proteo¬ 
glycan  content  may  also  be  exploited  to  formulate  an 
imaging  method  that  could  detect  only  a  small  loss  of 
proteoglycans  from  cartilage.  The  prolongation  of  the  Tx 
implies  that  Tx -weighted  images  would  indicate  a  greater 
degree  of  signal  loss  in  regions  of  cartilage  with  a  lower 
proteoglycan  content.  The  lengthening  of  the  T2s,  however, 
will  tend  to  increase  the  signal  from  degraded  cartilage. 
Short  echo-time  images  will  limit  the  increase  in  signal 
secondary  to  this  effect  while  still  exploiting  any  reduction 
in  r2f  that  will  act  to  decrease  the  signal.  Short  echo  times 
also  retain  the  adequate  signal  levels  critical  in  sodium 
imaging.  Single  quantum  sodium  images  of  cartilage  that 
has  undergone  serial  enzymatic  degradation  have  been 
presented  to  demonstrate  that  this  technique  may  detect 
the  loss  of  proteoglycans  from  the  cartilage.  This  indicates 
that  a  noninvasive  method  of  detecting  the  early  changes  of 
osteoarthritis  may  be  feasible. 

In  summary,  the  loss  of  proteoglycans  from  the  matrix  of 
articular  cartilage  may  be  measured  by  both  single  quan¬ 
tum  sodium  spectroscopy  and  imaging.  The  loss  of  even  a 
small  amount  of  proteoglycans  from  cartilage  also  leads  to 
a  prolongation  of  both  the  T2s  and  7\  of  the  sodium  NMR 
signal.  A  decrease  in  the  T2 f  is  also  detectable  with  greater 
amounts  of  proteoglycan  loss  from  the  cartilage.  These 
results  indicate  that  short-echo,  Ti-weighted  single  quan¬ 
tum  sodium  imaging  will  detect  the  early  stages  of  osteoar- 
thritic  change  in  cartilage.  Single  quantum  sodium  images 
of  cartilage  specimens  that  have  undergone  serial  proteo¬ 
lytic  enzymatic  degradation  were  presented  to  demon¬ 
strate  this  fact. 

REFERENCES 

1.  Freeman  MAR.  Adult  Articular  Cartilage,  2nd  ed.  England:  Pitman 
Medical  Publishing,  1979. 

2.  Lesperance  LM,  Gray  ML,  Burstein  D.  Determination  of  fixed  charge 
density  in  cartilage  using  nuclear  magnetic  resonance.  J  Orthop  Res 
1992;10:1-13. 

3.  Sapolsky  AI,  Kreiser  H,  Howell  DS,  Woessner  JF.  Metalloproteinases  of 


human  articular  cartilage  that  digest  cartilage  proteoglycan  at  neutral 
and  acid  pH.  J  Clin  Invest  1976;  58:1030-1041. 

4.  Woessner  JF,  Selzer  MG.  Two  latent  metalloproteinases  of  human 
articular  cartilage  that  digest  proteoglycan.  J  Biol  Chem  1984;  259:3633- 
3638. 

5.  Martel-Pelletier  J,  Pelletier  J-P,  Clouteir  JM,  Howell  DS,  Ghandur- 
Mnaymneh  L,  Woessner  JF.  Neutral  proteases  capable  of  proteoglycan 
digesting  activity  in  osteoarthritic  and  normal  human  articular  carti¬ 
lage.  Arthritis  Rheum  1984;27:305-312. 

6.  Ehrlich  MG.  Degradative  enzyme  systems  in  osteoarthritic  cartilage.  J 
Orthop  Res  1985:  3:170-184. 

7.  Ralph  Schumacher  H.  Ed.,  “Primer  on  the  Rheumatic  Diseases,”  10th 
ed,  Atlanta,  Ga;  Arthritis  Foundation,  1993. 

8.  Maroudas  A,  Evans  H.  A  study  of  ionic  equilibria  in  cartilage.  Connect 
TissRes  1972;  1:69-77. 

9.  Foy  BD,  Lesperance  LM,  Gray  ML,  Burstein  D.  NMR  parameters  of 
interstitial  sodium  in  cartilage,  in  Proc.,  SMRM,  8th  Annual  Meeting, 
Amsterdam,  1989.  p.  1108. 

10.  Bashir  A,  Gray  ML,  Burstein  D.  Sodium  Tj  and  T2  in  control  and 
degraded  cartilage:  implications  for  determination  of  tissue  proteogly¬ 
can  content,  “Proc.,  SMRM,  14th  Annual  Meeting,  Nice,  1995,”  p.  1896. 

11.  Jelicks  LA,  Paul  PK,  O’Byrne  E,  Gupta  RK.  Hydrogen- 1,  Sodium-23,  and 
Carbon- 13  MR  spectroscopy  of  cartilage  degradation  in  vitro.  JMRI 
1993;3:565-568. 

12.  Tyler  JA.  Cartilage  degradation.  Cartilage:  Molecular  Aspects,  {Hall  B, 
Newman  S,  editors),  Chapter  7,  CRC  Press,  Boca  Raton,  1991 . 

13.  Farndale  RW,  Sayers  CA,  Barnett  AJ.  A  direct  spectrophotometric 
microassay  for  sulfated  glycosaminoglycans  in  cartilage  cultures.  Conn 
TissRes  1982;9:247-248. 

14.  Jacard  G,  Wimperis  S,  Bodenhausen  G.  Multiple  quantum  NMR  spectros¬ 
copy  of  S  =  3/2  spins  in  isotropic  phase:  a  new  probe  for  multiexponen¬ 
tial  relaxation.  J  Chem  Phys  1986;  85(ll):6282-6293. 

15.  Bottomley  PA.  A  versatile  magnetic  field  gradient  control  system  for 
NMR  imaging.  J  Phys  E  Sci  Instrum  1981;  14:1081-1087. 

16.  Insko  EK,  Reddy  R,  Leigh  JS.  High  resolution,  short  echo  time  sodium 
imaging  of  articular  cartilage.  JMRI  1997;  7 (6):  1056- 1059. 

17.  Provost  TJ,  Hurst  GC.  Asymmetric  sampling  in  2D  FT  magnetic  reso¬ 
nance  imaging,  in  “Proc.,  SMRM,  5th  Annual  Meeting,  Montreal, 
1986,”  p.  265. 

18.  Macovski  A,  Noll  DC.  Homodyne  detection  for  improved  SNR  and 
partial  k-space  reconstruction,  in  “Proc.  SMRM,  7th  Annual  Meeting, 
San  Francisco,  1988,”  p.  815. 

19.  Oh  CH,  Hilal  SK,  Ra  JB,  and  Cho  ZH.  Faster  magnetic  resonance  imaging 
by  use  of  3-quarter  matrix  data,  Proc.,  SMRM,  6th  Annual  Meeting,  New 
York,  1987.  p.  455. 

20.  Stoyanova  R,  Kuesel  A,  Brown  T.  Application  of  principle-component 
analysis  for  spectral  quantitation.  J  Magn  Reson  1995;  Al  15:256-269. 

21.  Elliott  MA,  Walter  GA,  Swift  A,  Vandenbourne  K,  Schotland  JC,  Leigh 
JS.  Spectral  quantification  by  principle-component  analysis  using 
complex  single  value  decomposition,  Proc.,  1SMRM,  6th  Annual  Meet¬ 
ing,  Sidney,  1998,  p.  1855. 

22.  Bevington  PR.  Data  Reduction  and  Error  Analysis  for  the  Physical 
Sciences,  New  York:  McGraw-Hill;  1969.  pp.  237-239. 


